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For  system  performance  estimation  and  the  inter- 
pretation of  measured  data,  there  is  a need  to  predict  and 
simulate  the  statistical  properties  of  the  response  of  a 
horizontal ly-directive  acoustic  array  to  ambient  noise. 

In  current  use  in  the  Navy  there  are  a number  of  such 
"beam-noise"  models  which  predict  the  properties  of  noise 
caused  by  the  dominant,  prevailing  source  at  low  frequen- 
cies; surface  shipping.  The  objective  of  this  study  is 
to  review  these  models  in  some  detail  and  then  to  recommend, 
for  specific  applications,  approaches  which  utilize  the 
best  features  of  each  model. 

V 

Nine  Navy-sponsored  models  were  thus  investi- 
gated, and  all  were  found  to  be  based  on  the  fundamen- 
tal hypothesis  that  the  beam  noise  is  the  convolution  of 
the  array  beam  pattern  with  the  sum  of  intensities  from 
the  individual  ship  sources.  The  models  are  then  distin- 
guishable by  their  treatments  of  the  ship  sources,  the 
transmissibn  loss  (TL),  and  the  array's  response.  Two 
categories  wore  natural;  Analytic  models  which  calculate 
the  statistical  properties  of  noise  directly  from  those 
of  the  components  (source  level,  TL,  etc);  and  Brute- 
Force  models  which  use  simulation  or  Monte  Carlo  techniques. 
The  ability  to  calculate  statistics  over  all  possible  val- 
ues of  a parameter  (e.g.,  ship  locations  for  a given  den- 
sity, all  possible  source  levels),  or  over  values  which 
might  occur  in  a short  time  period  {aay  24  hours)  led  to 
a second  type  of  classification:  "grand"  and  "short-term" 
ensembles. 
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The  nine  models  have  been  reviewed  and  are  des 
cribed  in  terms  of  approach,  treatment  of  ships  and  TL, 
receiver  submodel,  implementation,  advantages,  and  short 
comings.  Tables  are  given  to  summarize  these  findings. 

For  various  applications  of  beam-noise  models, 
the  statistical  quantities  and  types  of  ensembling  re- 
quired have  been  proposed  (e.g.,  first-order  statistics 
with  full  ensembles,  or  multi-array  statistics),  and 
the  most  promising  models  and  techniques  are  identified 
for  each.  General  recommendations  for  a LRAPP  approach 
are: 


(a)  For  fully-ensembled  first-order  statistics, 
use  an  automated  version  of  the  USI  Analy- 
tic model. 

(b)  Use  the  BTL  Analytic  model  formulas  to 
obtain  bounds  on  temporal  statistics  and 
to  determine  dominant  parameters  for  the 
full  ensemble. 

(c)  For  details  of  fluctuations,  including 
effects  of  beam  patterns  and  Tt  variations, 
measurement  interpretation,  detection  stu- 
dies, use  the  best  of  the  Brute- Force  models 
or  a synthesis. 

(d)  For  multibeam  correlation,  use  BTL  formu- 
las for  bounds  and  a Bruce-Force  model 
for  details. 

(e)  Further  investitjation  is  required  in  a 
number  of  problem  areas  common  to  all  of 
the  models:  Ship  Information  (Inputs), 

Wind  Noise,  Model  Evaluation,  Statistics 
Appropriate  to  Dynamic  Array  and  Detector 
Simulation. 
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Section  1 


INTRODUCTION 


This  report  presents  a critical  review  of  a 
number  of  the  approaches  currently  used  within  the  Navy 
for  modeling  the  statistics  of  low-frequency  beam  noise. 

It  concentrates  on  the  underlying  assumptions,  the  type 
and  resolution  of  required  input,  and  the  extent  of  the 
statistics  modeled.  The  objective  is  to  arrive  at 
recommended  conceptual  models,  which  utilize  the  best 
features  of  those  reviewed,  for  specific  applications. 

’’Beam  noise"  in  this  paper  means  the  processed 
response  (output)  of  a "beamed  system"  to  ambient  sea 
noise.  A "beamed  system"  is  an  array  with  horizontal 
directivity  and  a bearoforroer  which  responds  to  plane  waves 
in  accordance  with  a "conventional"  lorizontal  beam  pattern 
(main  beams  and  sidelobes).  It  is  assumed  that  the  system's 
temporal  processing  includes  both  narrowband  filtering 
and  Incoherent  averaging,  so  that  the  system  output  is 
proportional  to  an  incoherent  average  of  noise  intensities 
in  a narrow  band.  Since  we  focus  on  low  acoustic  fre- 
quencies, the  models  reviewed  concentrate  on  the  fluctua- 
tions of  the  surface-ship  coi)iq)onent  of  noise. 
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Nearly  any  of  the  models  which  predict  mean  beam- 
noise  levels  (TASSRAP,  RANDI,  FANM,  etc.)  can  also  in  theory 
be  used  to  predict  statistics  or  time  histories  via  multiple 
runs  with  inputs  varied  appropriately  (ship  densities  or 
locations  or  source  levels,  transmission  loss,  beam  pat- 
terns). However,  this  review  considers  only  those  models 
which  are  designed  and  structured  to  yield  the  statistical 
properties  of  noise: 

A-  USI  Model  (Underwater  Systems,  Inc.) 

B-  BTL  Model  (Bell  Laboratories) 

C-  BBN  Model  (Bolt  Beranek  and  Newman,  Inc.) 

D-  WAGNER  Model  (Wagner  Associates) 

E-  NABTAM  Model  (ORI,  USI,  NORDA,  et  al) 

F-  DSBN  Model  (Science  Applications,  Inc.) 

G-  BEAMPL  (NORDA) 

H-  SIAM  I (NRL) 

I-  SIAM  II  (NRL) 

Each  of  these  is  described  in  Section  2,  while  Section  3 
summarizes  their  attributes  and  makes  recommendations  for 
the  synthesis  of  a model  for  LRAPP  applications. 

Since  many  of  the  models  are  without  up-to-date 
documentation,  the  author  (or  party  responsible)  was  asked  to 
review  a draft  copy  of  the  section  of  this  report  which 
dealt  with  his  particular  mode! . The  reviewers  were  as 
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follows: 


Model 


Reviewers 


R.  L.  Jennette 
Underwater  Systems,  Inc. 


J.  Goldman 
Bell  Laboratories 


U.  Moll 

Bolt  Beranek  and  Newman,  Inc. 


B.  J.  McCabe 

Daniel  H.  Wagner,  Associates 


J.  J.  Cornyn 
NORDA 

and 

E . J . Moses 

Operations  Research,  Inc. 


C.  W.  Spofford 

R.  G.  Sticglitz 

Science  Applications,  Inc. 


S.  C.  Wales 

Naval  Research  Laboratory 
and 

S.  W.  Marshall 
NORDA 
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Written  comments  were  received  from  each  reviewer, 
and  corrections  have  been  incorporated  in  the  text  where 
appropriate.  The  authors  of  this  report  appreciate  this 
cooperative  effort  to  make  the  model  descriptions  a.s 
accurate  as  possible. 

The  remainder  of  this  introductory  section  des- 
cribes in  brief  the  various  fluctuation  mechanisms,  defines 
beam-noise  fluctuation  parameters  according  to  time  scales 
and  ensembling,  and  identifies  the  general  approaches  to 
n^deling  noise  statistics. 


1.1  FLUCTUATION  MECHANISMS.  TIME  SCALES  AND  ENSEMBLES 

Beam  noise  fluctuations  are  the  result  of  any 
of  a number  of  mechanisms.  At  low  frequencies  the  dominant 
prevailing  sources  of  ambient  noise  are  surface  shipping 
and  wind  action  on  the  sea  surface,  so  that  the  mechanisms 
include: 

(A)  surface  ship  movements  on  and  off  the  main 
beams,  and  in  the  sldelobes. 

(B)  variations  in  average  transmission  loss 
(TL)  from  ship  source  to  receiver  causr^ 
oy  source/recelver  (S/R)  motion  in  a 
variable  environment  and  changes  in  S/R 


separat ion , 


(C)  fluctuations  in  the  detailed  TL  caused  by 
environiaental  variability  and,  for  narrow- 
band  noise,  S/R  movement  through  the  multi- 
path  interference  field, 

(D)  array  response  variation,  from  array  dis- 
tortion or  signal  processing  artifacts. 

(E)  surface-ship  source  level  fluctuations, 
both  from  changes  in  ship  aspect  and  from 
fluctuations  in  the  short-time  averaged 
radiated  noise  itself, 

(P)  variations  in  wind-generated  noise,  caused 
by  changes  in  the  environment , wind  speed 
and  direction. 

None  of  the  models  reviewed  explicitly  treats  all  of 
these.  If  we  disregard  time  scales,  system  details,  and 
ensembles  for  the  moment,  then  the  mechanisms  believed 
to  be  most  io^ortant  and  modeled  most  often  are  ranked  in 
order  above.  The  fluctuations  and  directionality  of  wind 
noise  <P)  may  be  very  important  in  some  circumstances, 
but  little  is  known  about  them  and  the  mode's  considered 


do  not  predict  them. 
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In  order  to  properly  compare  and  assess  the  beam 
noise  models,  we  must  better  define  the  quantities  to  be 
estimated  - the  beam-noise  statistics.  Assuming  that  the 
model  output  is  to  be  used  to  assess  system  detection  per- 
formance and  avoiding  a treatise  on  detectors,  we  simply 
state  some  typical  applications  of  beam-noise  description 

- instantaneous  detection  probability 

- cumulative  detection  probability  in  several 
hours  or  a day 

- distribution  of  holding  times  or  non-holding 
times  in  a day. 

and  some  corresponding  beam-noise  "statistics”: 

- distribution  of  beam-noise  level  over  a short 
time  period  (minutes  to  an  hour) 

* 

- statistical  dependence  of  noise  at  points  in 
time  separated  by  minutes  to  hours  or  a day. 

- distribution  of  beam  free  times,  over  hours 
or  a day. 

- distribution  of  time  intervals  for  which  noise 
is  below  a threshold,  over  an  hour  or  a day. 


owinp^tatr 


Now,  each  of  these  statistical  distx'ibutions  or  correlation 
functions  is  based  on  a temporal  "average"  over  the  time 
period  specified.  For  example,  if  beam  noise  was  observed 
as  a time  series  N(t),  then  its  distribution  function  over 
period  < t < t2  is  simply  the  distribution  function  of 
the  set  of  points  |N(t)|  < t < tgf  with  each  point  in 

time  (t- , t„)  equally  likely;  F„(X)  = PtN(t)  < X for 
ti  < t < tg] . In  estimating  such  a distribution  function 
with  a noise  model,  the  environmental  and  noise-source 
conditions  can  change  only  as  much  as  would  be  observed 
in  that  time  period  (say,  tj^  t < tg).  In  most  applications, 
the  system  performance  is  to  be  described  for  an  ensemble 
of  conditions,  e.g. , over  all  ship  distributions  consis- 
tent with  some  average  densities,  or  over  transmission 
loss  conditions  occiirring  in  a week  or  throughout  a geo- 
graphic area.  The  appropriate  quantities  to  be  estimated 
by  the  noise  model  are  then  the  ensemble  of  the  time  sta- 
tistics discussed  above,  e.g.,  the  set  of  distribution 
functions  {Fj^  (^)|l  generated  when  the  ensemble  of  condi- 
tions is  considered  (denoted  by  an  index  set  I).  If  short- 
term statistics  are  desired,  then  it  is  usually  a mistake 
to  "average"  over  the  ensemble.  Again,  as  an  example, 
suppose  that  over  given  conditions  i,  the  noise  distribu- 
tion function  F„  (X)  is  normal  with  mean  y.  and  variance 


ii>iW^lfi>lMl(i.iiliiM('iMi«iiiUi^*^ 
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0^  . Those  are  the  quantities  which  determine  detection 
probability.  If  the  ensemble  of  conditions  is  indexed  by 
a finite  set  i«I,  then  the  ensemble  distribution  function 
Fj^(X)  is  found  as 

P{^(X)  *»  PlN^(t)<X  over  all  t and  i], 
n 

- 2;  PIN,  (t)<XU-ijjl.P(l=i^) 

1^-1  o 

“ ^ E PIN.  (t)<X] 

io-1  " ' 

i.e.,  , is  the  average  of  the  distribution  functions 

of  the  , is  no  longer  normal , and  may  have  a very  large 

2 ^ 2 
0 - larger  than  any  of  the  individual  0^  . This  type  of 

result  can  be  very  misleading  in  the  calculation  of  detec- 
tion probabilities. 

The  usual  px'obabilistic  treatment  of  the  time- 
average/ensemble-average  problem  is  to  use  assumptions 
about  the  ergodicity  (and  hence  stationarity ) of  a noise 
process  N(t,i)  (a  time  aeries  for  each  i€l,  a random 
variable  over  the  ensemble  i<I  for  each  fixed  t).  In  that 
case,  "long"  time  averages,  over  many  independent  samples, 
and  ensemble  averages  are  equivalent: 


’fpnnsurtfKr*- 


There  may  be  justification  for  the  ergodicity  assumption 
for  noise  over  some  specific  ensembles,  but  in  general 
for  the  problems  of  interest  here  there  is  no  basis  for 
it  - and  even  if  there  were,  the  time  average  would  be 
over  very  long  periods. 


To  conclude  we  list  below  the  typical  time  scales 
for  the  fluctuation  mechanisms  identified  at  the  beginning 
of  the  subsection: 


Mechanism 


Typical  Time  for  a 
Significant  Change 


• Ships  on  and  off 
beams 


Minutes  for  Nearby  Ships , 
Hours  to  Days  for  Distant 
Ships 


• Average  TL  variations  Hours 


• Detailed  TL  fluctua-  Minutes  to  Hours 
tions 


• Array  response  varia-  Minutes  to  Hours 
tions 


• Ship  source  level  Hours  to  Days 

variations 


• Wind  noise  variations  Hours  to  Days 
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As  mentioned  earlier,  all  of  the  models  reviewed 
here  concentrate  on  the  ambient  noise  component  resulting 
from  surface-ship  sources  and  treat  wind  noise  only  as 
an  additive,  empirical  terra.  Hence,  those  attributes 
which  distinguish  one  model  from  another  are  the  input 
requirements,  the  model  output,  the  methods  of  calculating 
the  noise  properties,  and  treatments  of  ship  sources,  TL 
and  array  response. 


where  J(t)  is  the  number  of  ships  at  time  t, 

!SLj(t)  is  the  .^ource  intensity  of  the  j-th  ship 
> at  time  t, 

Tj(t)  is  the  transmission  ratio  for  ship  j at 
time  t (i.e.  TL  » -lOlog(Tj)), 

AGj(t)  is  the  arr^y  response  for  the  arrivals 
from  ship  j ."t  time  t. 

♦NABTAM  is  an  exception.  It  calculates  a directional, 
time-independent  wind-noiso  component. 


When  all  of  these  quantities  are  known,  deterministic  func- 
tions of  time,  then  N(t)  is  directly  available  from  (1-1); 
this  is  the  way  that  most  average-noise  models  work,  (i.e., 
"point"  models  such  as  TASSRAP  or  RANDI). 

The  form  of  (1-1)  is  not  meant  to  imply  that  the 
individual  factors  are  necessarily  uncoupled.  In  fact, 
in  some  of  the  models  the  source  intensity  and  array  re- 
sponse depend  on  the  transmission  ratio  (e.g.,  multipath 
angles  for  a geometric  treatment).  In  the  same  vein,  a 
precise  model  of  the  temporal  behavior  of  SL  would  have  to 
reflect  the  transmission  travel  times.  Finally,  for  some 
applications  it  might  be  necessary  to  perform  the  summation 
of  intensities  in  (1-1)  on  a phased  basis. 

Consider  now  two  basically  different  approaches 
to  computing  beam-noise  statistics. 

Analytic  Approach;  An  Analytic  approach  is  one  in  which 
the  variables  of  (1-1)  are  treated  as  random  processes 
with  well-specified  statistical  properties,  such  as  k- 
diraensional  distribution  functions.  Standard  probabilistic 
techniques  (convolution,  characteristic  functions,  etc.) 
are  then  used  to  calculate  statistical  properties  of  N(t) 
from  those  of  SL,  T,  AG  and  J.  For  example,  if  all  the 
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variables  in  the  sum  of  (1-1)  are  independent,  then  the 
mean  value  of  the  noise  intensity  can  be  found  directly 
from  the  mean  values  of  the  terms  in  the  sum.  Such  an 
"analytic"  approach  is  attractive  sinco  it  can  produce  an 
efficient  model  for  predicting  selected  noise  statistics 
and  for  identifying  the  mechanisms  which  are  important. 

On  the  other  hand,  the  calculations  can  become  complicated 
when  higher-order  moments  or  distribution  functions  or 
"short-term"  statistics  as  described  above  for  N(t)  are 
required.  Moreover,  some  of  the  assumptions  required  to 
make  the  calculations  tractable  may  be  unrealistic,  and 
not  all  of  the  statistical  properties  of  N may  be  available. 
It  is  the  clever  implementation  of  the  approach  which  makes 
it  useful. 


The  BBN,  WAGNER,  BTL,  and  USI  models  employ  an 
Analytic  approach. 

Brute-Force  Approach:  The  second  approach  can  be  viewed 
as  a special  case  of  the  Analytic  approach,  in  which  the 
statistical  properties  of  N(t)  are  evaluated  numerically  - 
by  calculating  realizations  of  N(t)  for  specific  values  of 
SL,  T,  AG  and  J.  All  of  the  quantities  of  (1-1)  are 
treated  as  known  and  deterministic  for  a single  realiza- 
tion of  N(t)  over  some  time  period,  say  0 ^ t < T.  In 
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particular,  ship  tracks  and  source  levels,  transmission 
loss  versus  range,  angles  and  time,  and  the  array  response 
are  all  specified  over  tO,TJ . This  might  be  the  case  when 
the  results  of  a measurement  exercise  are  to  be  simulated 
with  the  model.  If  there  is  to  be  ensembling  over  any 
of  these  variables,  then  additional  realizations  of  N(t) 
are  generated  from  new  ship  tracks  or  TL  or  whatever. 

This  is  usually  called  a Monte  Carlo  simulation  in  that 
the  description  of  N(t)  over  the  ensemble  of  conditions 
is  found.  It  is  not  (necessarily)  a Monte  Carlo  method 
in  the  ordinary  statistical  sense  that  the  individual 
estimates  are  combined  and  the  combination  converges  to 
some  ensemble  descriptor  (e.g.,  the  distribution  of  the 
mean  value). 

The  Brute-Force  approach  can  use  all  of  the 
information  about  ships,  TL,  and  array  response  available 
and  can  yield  just  about  any  noise  statistic  including  the 
"short-term"  properties  mentioned  above.  On  the  other 
hand,  it  can  be  very  time-consuming,  often  requiring 
many  realizations  to  cover  the  ensemble  conditions  and 
roust  rely  on  special  analysis  routines  to  summarize  or 
process  or  Interpret  the  large  amount  of  output. 
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BEAUPL  (NORDA).  DSBN  (SAI).  SIAM  I and  II  (NRL). 
and  NABTAH  (NORDA)  are  Brute-Force  models. 

In  summary,  each  of  the  two  approaches  has  ad- 
vantages and  disadvantages,  vrhose  Importance  depends  on 
the  application  of  the  model  predictions.  The  term 
"Brute-Force"  is  not  meant  to  suggest  that  such  models 
are  unsophisticated  or  lack  mathematical  insight.  Nor 
should  the  term  "Analytic"  convey  more  than  the  notion  that 
the  calculation  is  performed  without  direct  realizations 
of  the  random  vax'iable  (a  numerical  calculation  of  a 
characteristic  function  and  ics  inversion  can  be  as  tedious 
as  a Monte-Carlo  computation). 
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Section  2 


REVIEW  OF  BEAM- NOISE  MODELS 


This  section  describes  each  of  the  nine  beam- 
noise-statistics  models,  concentrating  on:  the  input  re- 
quirements, treatment  of  ships,  TL,  array  response,  com- 
puter calculations,  output  quantities  and  analysis,  status 
of  computer  codes  and  documentation.  The  order  in  which 
the  models  are  presented  and  the  number  of  pages  devoted 
to  a single  approach  were  convenient  to  the  author  and 
have  no  other  significance. 

An  attempt  has  been  made  to  utilize  consistent 
notation  for  the  model  descriptions: 


N or  I : 

noise  intensity 

SL: 

source  intensity 

T: 

transmission  ratio 

TL: 

transmission  loss  (10  log  T) 

AG: 

array  power  response 

I 


R:  range  from  array 

4)  or  li*:  bearing  from  array  or  angle  from  array 
broadside* 

6:  vertical  angle  of  arrival 

J or  K:  number  of  ships 

X:  mean  of  Poisson  variable 

V : speed 

u:  mean 

o:  standard  deviation 

f;  density  function 

characteristic  function* 

P:  probability 

E:  expected  value 


Var(*)-  variance 

Finally,  a few  definitions  are  reviewed. 

For  random  variables  X, . X , the  multi 

X ^ n 

variate  density  function 


r 

I 

[ 

I 


*There  should  be  no  confusion  in  using  for  the  charac- 
teristic function  and  for  bearing  angle. 
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Mrill  be  called,  according  to  standard  terminology,  an 
n-dimensional  or  n-tfa  order  density  function.  A set  of 
random  variables  indexed  by  t,  is  a random  process 

or  stochastic  process.  Its  n-dimensional  density  is  the 
multivariate  density  of 

X(tj).  XCtg).  ....  X(t^) 

for  indices  tj^,  tg.  ....  t^.  A sample  path  or  realiza- 
tion from  X(t)  consists  of  a sample  from  X(t)  for  each 
index  t.  An  n-tfa  order  statistic  for  X(t)  is  a statistic 
which  depends  (nontrlvlally ) on  the  n-dimenslonal  density. 
Hence,  the  mean,  variance,  skewness,  median,  and  one- 
dimensional  density  for  X(t)  are  first-order  statistics, 
while  the  autocorrelation  function  of  X(t)  is  a second- 
order  statistic. 

The  characteristic  function  for  a random 
variable  X is 

♦(«)  ■ E(oxp( ib>X)) . 
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USI  MODEL 


2.1.1 


Background 


Name:  (USI  Array  Noise  Model.  Version  1)  (USI) 


Developer:  Underwater  Systems,  Inc., 

R.  L.  Jeanette  and  E.  L.  Sander 


Sponsor:  Applied  Physics  Laboratory. 
Johns  Hopkins  University 


Previous  Applications:  APL/JHU  and  LRAPP  studies 


Published  Documentation:  Reference  A-5  , A-6 


2.1,2 


General 


The  USI  model  uses  an  Analytic  apni’oach.  It 
numerically  estinuxtes  the  ensemble  and  time-averaged,  one- 
dimensional statistical  distribution  lunction  of  beam 
noise  using; 


Time-averaged  ship  densities  for  %'ach  of 
a set  of  contiguous  get>graphical  regions 


• P<.>isson  m(.>del  for  ship  counts  in  each 
geographic  region 


• Uniform  naidel  for  ship  UK'ations  In  a 
region 
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• Input  averaged  TL  versus  range 

• Input  TL  fluctuation  distribution 

• Input  ship  source  levels  and  fluctuation 
distributions 

• Input  asimuthal  array  beam  pattern 

• Incohei’t^nt  addition  of  ship  contributions 

Preparation  of  the  data  and  the  formula  (1-1)  lead  to  a 
random- variable  model  for  noise  intensity  (1)  at  the  array 
output  of  the  form 


(A-l) 


where  '.I  (nli  are  random  variables  reprt>senting  the  possL 
* k ' 

ble  intensity  contributions  from  single  ships,  are 

Poisson  variables  representing  the  corresponding  ship 
counts,  and  K Is  the  number  of  distinct,  per-ship,  inten- 
sity distribut  ions.  the  assumi>tion  is  that  {lj^(n}(and 
{ { form  a set  of  independent  variables,  and  that  Jjj(n) 
has  the  same  distribution  for  all  n when  k is  fixed.  Ute 


computer  code  numerically  estimates  the  distribution  func- 
tion (with  1-1/2  or  3 d0  resolution)  from  equation  <A-1) 
using  Input  mean  values  of  {Jj^j  and  consistently  precise 
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A typical  computer  run  (execution  only)  costs 
less  than  a dollar,  and  the  code  produces  supplemental 
information  about  which  ships  are  dominant  in  determining 
the  mean  and  variance  of  the  noise.  Only  first  order 
statistics  are  available.  The  USI  model  differs  from  the 
others  in  its  fast  numerical  approach  to  estimate  the 
ensembled  one'dimensional  distribution  function  for  beam 
noise. 


2.1.3 


Model  For  Ships 


■ I 


At  present  the  ship  fields  and  intensity  esti- 
mates are  constructed  manually.  The  usual  ship  data  are 
average  ship  densities  for  ocean  regions,  say  x l^ 
cells  or  5^  x cells.  Itthin  each  region  the  actual 
number  of  ships  is  treated  as  a Poisson  variable,  inde- 
pendent from  one  region  to  another  (at  a fixed  time), 
and  the  ships  are  assumed  to  be  distributed  unifoiinly  in 
the  region. 


The  IHiisson  assumption  is  consistent  with 
classical  probability  models  for  distributing  points  in 
space  when  the  points  are  not  allowed  to  cluster  (see, 
e.g..  Feller  (1957).  p.  149).  That  assumption  is  also 
critical  to  the  formulation  of  the  USI  approach,  which  uses 
a special  property  of  the  Poisson  variable;  if  the  ship 
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count  in  an  area  is  a Poisson  variable  with  mean  M and  the 
area  is  subdivided  into  smaller,  non-overlapping  subregions, 
then  the  shin  counts  in  the  subregions  are  independent 
Poisson  variables  with  means  proportional  to  the  subregion 
areas  and  with  the  sum  of  the  means  equal  to  M.  In  the 
USI  approach,  regions  are  broken  up  and  regrouped  accord- 
ing to  the  contribution  of  a ship  in  the  region  to  the 
noise  intensity.  The  Poisson  assumption  allows  the  model 
to  treat  the  ship  count  for  the  grouped  regions  as  a 
Poisson  variable  with  the  logical  mean  value.  Note  here 
that  there  is  no  explicit  temporal  variation  of  ship  loca- 
tions or  properties. 


Ships  in  a region  are  assigned  a source-level 
density  function,  with  3 dB  (or  1-1/2  dB,  if  desired) 
resolution.  Construction  of  these  functions  is  part  of 
the  manual  preparation  of  the  input. 


The  final  step  in  preparation  of  the  ship  data 
is  to  group  ships  according  to  the  distribution  of  inten- 
sity at  the  receiver  output.  This  is  accomplished  by 
incorporating  the  distribution  function  for  TL  (and  even 
array  response)  ir  the  source  level  distribution  as  a 
function  of  the  location  of  the  region,  both  in  range 
and  azmimuth.  The  result  then  is  that  for  a given 
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azimuthal  sector,  there  is  a set  of  intensity  variables 
|lj^}  each  with  distribution  derived  as  the  sum  of  inde- 
pendent variables  for  source  level,  transmission  loss, 
and  array  response,  and  there  is  a corresponding  set 
of  Poisson  ship-count  variable  |jj^[  such  that  Jj  repre- 
sents the  number  of  ships  in  the  sector  with  intensity 


2.1.4 


Model  for  Transmission  Loss 


Transmission  loss  and  fluctuation  distribution 
are  inputs  to  the  USI  model,  and  are  in  fact  used  only 
in  the  input  preparation  described  above.  To  be  precise, 
the  transmission  ratio  (intensity  loss)  is  treated  as  a 
deterministic,  smooth  or  averaged  function  of  range  to 
which  is  added  a fluctuation  term.  The  smooth  transmission 
is  used  to  group  the  ship  regions  according  to  intensity 
contribution  to  the  noise,  while  the  fluctuation  term  is 
added  as  an  independent  variable  to  the  appropriate  ship 
source  intensity  fluctuations. 


When  the  resolution  of  the  source  level  and 
output  intensity  distribution  is  X dB  (3  or  1-1/2),  then 
the  average  values  and  distributions  of  the  transmission 
ratio  are  normally  estimated  with  similar  precision. 


r 

4i<- 

i 

i 

I 


the  entire  population  of  ship  distributions  (for  the 
given  average  density  field).  The  corresponding  time 
period  in  replication  of  ship  locutions  required  to 
realize  this  amount  of  variation  is  probably  on  the 
order  of  days  or  longer.  The  need  (and  ability)  to 
model  the  details  of  short  term  fluctuations  at  the 
processor  seems  then  to  be  diminished.  The  underlying 
assumption  made  in  the  USI  approach  is  that  "the  noise 
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process  has  a long  correlation  time."  More  on  this 
subject  can  be  found  below. 

2.1.6  Details  of  the  Calculation 

Most  of  the  previous  discussion  concerns  the 
manual  preparation  of  the  input  to  the  USI  model.  The 
actual  computer  routine  operates  on  the  input  data: 

- intensity  for  group  k,  a random  variable 
with  given  distribution  function, 

“ number  of  ships  in  group  k,  a Poisson  variable, 
to  estimate  the  distribution  function  for  the 
received  noise  intensity: 

K 

I - 2 ( £ • (A-1) 

k-1  n»l 


There  are  several  ways  in  which  this  distribution  could 
be  found,  and  it  is  USI's  approach  which  makes  the  model 
unique.  Consider  then  the  "classical"  approach  as 
contrasted  to  USI's. 
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2. 1.6.1  CLASSICAL  APPROACH 

Consider  first  the  properties  of  I and  the 
"classical"  way  of  estimating  its  distribution  function. 
Suppose  that  the  Poisson  variables  Nj^  have  mean  Then 

e'^k(X  )*" 

p(j  » m)  = ^ , (A-2) 

* ml 

E(Jjj)  = Var  (J^)  = X^.  (A-3) 

Let  Fj^(x)  denote  the  distribution  function  for  ■=  Ijj(n) 
and  let 

<l>jj(w)  = Etexp(iwlj^)]  (A-4) 

denote  the  characteristic  function  for  Ij^. 


Now,  the  inner  summand  of  (A-1), 


Jk 


n“l 


(A-5) 


can  be  viewed  as  a fixed  time  point  of  a compound  Poisson 
process  (see  any  book  on  stochastic  processes  such  as 
Papoulis  or  Parzen).  The  independence  of  the  Ijj(n)  and 
Jj^  allow  for  the  immediate  calculation  of  moments  of  Sj^: 
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E(SJ  ) “ E = n)  . P(Jj^  » n) 

n=o 

00 

= J])  EClf).  n*P(J.  =n) 
n=0  * 

» E(I  J) -ECJj^).  (A-6) 

In  particular, 

E(Sj^)  = Xj^E(Ij^) 

E(Sj^2)  = = Var(Sj^)  + tE(Sj^)l^  . (A-7) 


Furthermore,  the  characteristic  function  for  Sj^  is  given 
by 


^ (w)  = E(exp(iw  E 

k k^l  * 


= E E(iwl.  •n)-P(J.  =n) 
n=0  * “ 


E )e~^k^^k^^^"^\  (A-8) 


n»0 


HT“ 


The  independence  of  the  and  |jj^|  imply 

that  of  the  jSj^l  so  that  the  properties  of  I can  be  found 
from  (A-6)  and  (A-8): 
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- K „ K „ 

m ) » ZE(Sn  - EE(IJ)E(N.) 
k-1  * k»l  * 

(A-9) 

K 

Ed)  = 

k«l  * 

(A-10) 

2 K n 

E(i  ) “ E 

kS^l  ^ * 

(A-11) 

K 

<|>(w)  * E(exp(iwl)]  = n E(exp  iwS.  ) 

k=l  * 


n (w) 

k=l  ®k 


k»l 


exp( 


(A-12) 


(A-13) 


Hence,  all’  of  the  moments  of  I can  be  calculated  directly 
from  those  of  and  while  the  characteristic  function 
of  I is  given  as  a simple  function  of  that  of  in  (A-12) 
or  (A-13).  The  density  function  for  I is  then  the  inverse 
Fourier  transform  of  4), 


dw. 


f(x)  « ♦(w)e"^'*'^ 


(A-14) 
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Formulas  (A-13)  and  (A-14)  combine  to  yield  an 
algorithm  to  compute  f.  In  fact,  4>jj(w)  could  be  calculated 
with  an  FFT  (or,  if  were  Gaussian,  in  closed  form),  as 
could  the  inverse  transform  of  (A-14).  When  great  accuracy 
is  not  required  and  the  distributions  are  assumed  to  be 
smooth,  32  or  64  point  FFT's  could  be  used  with  time  of 

_3 

about  3 X 10  seconds  each. 


it  ^ 


A related  approach  to  finding  the  properties 
of  I is  to  evaluate  its  moments  directly  from  (A-9)  and 
then  to  use  an  Edgeworth  or  other  moment-expansion. 


The  developer  of  the  USI  model  has  found  an 
alternative  to  these  which  is  much  more  efficient  and 
which  readily  yields  information  about  the  dominant 
factors  affecting  the  noise  statistics.  It  is  outlined 
below. 


2. 1.6. 2 USI  Approach 


The  author  of  the  USI  model  has  described  the 
approach  as  follows  (details  appear  in  Reference  A-5). 


The  intensity  variable.  I,  is  viewed  as  in 


2-14 


(1-1): 
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I = E SL^*T.*AG.  = E W 

ships  ^ ^ ships  ^ 


(A-15) 


with  SL^,  and  AG^  random  variables.  The  grouping  of 
ships  in  domains  Dj^  such  that  a single  ship  in  the  k-th 
domain  has  "weight”  given  by  the  random  variable  Wj^  yields 


SI 


K Jk 

I = E ( E *k>  . 

k*l  j-1 

where  is  the  number  of  ships  in  Wj^  and  K is  the  number 
of  domains.  Now  is  assumed  to  be  a Poisson  variable 
and  the  author  writes 


I » E 


k "k  * 


(A-IG) 


This  intensity  is  then  treated  as  "the  sum  of  weighted 
Poisson  variables." 


To  numerically  estimate  the  distribution  func- 
tion for  I of  (A-16)  the  distribution  functions  for  the 
and  for  I are  discretized  in  3 dB  (or  1.5  dB)  steps 
and  a recursive  formula  is  used  to  determine  the  impor- 
tant contribution  to  I,  from  the  lowest  intensities,  in 
order,  to  the  highest.  This  procedure  allows  the  con- 
volution to  be  calculated  by  a simple  recursion,  and 
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limits  execution  costs  to  less  than  one  dollar  per  com- 
puter run.  The  number  of  3-dB  distribution  "bins”  can 
be  50  or  more,  providing  a range  of  noise  levels  of 
over  150  dB. 

2.1.7  Output  and  Analysis 

The  output  of  the  USI  model  consists  of 

• the  distribution  function  for  the  noise 
intensity  I,  v?ith  3 dB  (or  1.5  dB)  reso- 
lution. 

• an  analysis  of  how  each  geographical  region 
contributes  to  the  distribution,  including 
the  mean  and  variance. 

The  value  of  the  latter  is  to  help  to  focus  on  the  impor- 
tant ships  so  that  detailed  fluctuation  analysis  can  be 
performed.  We  note  again  that  the  distribution  function 
represents  a time  and  ensemble  average  over  all  ship  loca- 
tions consistent  with  the  input  densities  and  Poisson 
assumption . 


I 

I 

r 


7' 


I 

I 


% 

i 
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2.1.8  Computer  Implementation 

The  USI  program  is  coded  in  ANSI  Standard 
FORTRAN,  and  is  presently  executed  on  CDC  6000-series 
machines  and  the  Varian  620/L  minicomputer.  Core  re- 
quirement is  20K,  and  running  time  is  about  one  second 
or  less  on  the  CDC  6000. 

2.1.9  Evaluation 

Output  of  the  USI  model  has  been  compared  with 
the  noise  distribution  data  of  Reference  A-1  and  showed  favor- 
able agreement,  although  there  are  questions  in  the  inter- 
pretation of  the  data. 


2 . 1 . ^ 0 Significant  Advantages  and  Disadvantages 

1 

i 

I 


The  primary  drawback  for  operation  of  the  USI 

i 

, model  is  in  the  manual  preparation  of  the  shlp/intenslty 

groups.  Ship  cknsities  must  be  converted  to  azimuthal- 

:| 

* sector /range  bins  in  order  to  group  regions  by  Intensity 

r and  beam: 


i 
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USI  is  currently  Investigating  the  potential  for  utilizing 
FANIN,  a computer  routine  used  to  prepare  precisely  such 
ship  information  for  the  FANM  noise  model  at  NORDA  320. 

Other  shortcomings  include  the  fact  that  the  USI 
model  does  not  predict  ti^mporal  statistics.  Also,  the  key 
assumption  of  Poisson  ship  counts  has  not  been  thoroughly 
tested.  The  details  of  the  calculation  of  the  density 
function  have  not  as  yet  been  released,  so  that  a critical 
review  of  the  numerical  routine  (based  on  Equation  (A-16)) 
is  impossible  at  this  time.  Finally,  we  note  that  this 
model,  and  most  of  the  analytic  models,  predict  the  sta- 
tistics of  noise  over  an  ensemble  corresponding  to  ship 
movements  over  a long  time  period  (days).  Such  a predic- 
tion can  be  of  great  value  for  certain  applications,  but 
requires  careful  interpretation. 

On  the  positive  side,  we  note  that  the  USI  com- 
puter routine  is  oytreiaoly  fast  and  inexpensive  to  run. 
orders  of  magnitude  faster  than  the  other  models  reviewed 
here.  As  a production  tool,  with  automated  ship/Tt  input 
and  processing  of  geographic  domains,  the  approach  could 
be  extremely  valuable  for  a number  of  applications. 
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BTL  MODEL 

1 

2.2.1 

\ 

Background 

Name : ( BTL  Noise 

1 

Developer:  Bell  ! 

I 

Sponsor:  Naval  12 

Cod®  o20  and  PME-124 

Previous  Applications:  Analysis  of  surveillance- 

system  noise  data 

Published  Documentation:  References  B-1 , 8-2. 
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The  BTL  model  employs  an  analytic  approach.  It 
was  designed  for  the  surveillance  case  and  hence  to  pro- 
vide a statistical  description  of  shipping  noise  for  a 
(horizontal)  beamed  system  at  fxequencies  in  the  range 
25-150  Hz.  The  basis  of  this  and  the  other  noise  models 
IS  equation  (1-1),  but  the  key  to  the  BTL  approach  is 
the  assumption  that  ships  arrive  randomly,  according 
to  a Poisson  rule.  In  that  case  the  noise  intensity  can 
be  represented  as  a generalized  shot  noise  pi'‘ocess  (see. 
G.g.,  Parzen  (1862)).  Completely  characterized  by: 


i 

I 


i 


r 

r 

f 


I 
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• shipping  lanes* 

• mean  ship  arrival  time  on  each  lane 

• probability  distribution  of  ship  velocity 
and  source  level 

• probability  distribution  of  range/orienta- 
tion of  lanes  relative  to  the  beam  (or 
equivalently,  CPA  range  and  orientation) 

** 

• the  deterministic  transmission  loss,  array 
response , and  geometry 

Although  the  model  can  treat  arbitrary  shipping  scenarios, 
the  shot-process  formulation  has  been  exploited  in  special 
cases  to  yield  simple  analytic  expressions  for  the  statis- 
tics and  an  easy  identification  of  critical  factors 


affecting  them. 


Ctoldman  has  shown  that  the  noise  level  dis- 
tributions (in  dB's)  of  all  dimensions  tend  asymptotically 
(for  large  numbers  of  ships  or  for  "slow  ships  or  for  dis- 

« . 

tant  ships)  to  joint  normal  distributions,  then  such  a 
■r-  situation  is  valid,  the  approach  leads  to  an  especially 

r simple  analytic  model.  When  this  is  not  the  case,  calcu- 

lation of  density  functions  Is  performed  via  Fourier 

♦ ^ 

1 inversion  of  characteristic  functions. 

t *A.»i  will  be  discussed  below,  the  characterization  of  ship 

1 traffic  evolutions  by  lanes  is  necessary  only  to  simplify 

cuiculat ions. 

r •• 

f See  Reference  B-3. 
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Under  similar  assumptions  about  Poisson  ship 
arrivals,  the  BTL  approach  has  been  extended  to  the  multi- 
beam and  multi-array  cases  (Ref.  B-2).  Of  special  interest 
are  the  joint  density  functions  and  crosscorrelations  of 
noise  on  different  beams.  As  before,  the  ships  drive 
the  results  (TL,  SL,  AG  are  deterministic  and  time-inde- 
pendent), analytic  expressions  have  been  obtained  for 
special  cases,  and  the  joint  densities  are  shown  to  asyrap- 
totJcally  approach  normal  densities  (in  dB's). 

The  amount  of  input  preparation  and  computer 
expense  depend  greatly  on  the  particular  ship  scenario 
(e.g.,  number  of  lanes),  the  details  of  the  TL,  array 
response  and  source  functions,  and  the  type  of  output 
statistics  desired  (mor  ?nts,  1-diraensional  density, 
autocovariance,  multidimensional  density). 


* 2.2.3  Model  for  Ships 

r 

The  general  formulation  of  the  BTL  model  can 
f treat  an  arbitrary  collection  of  random  (straight)  ship 

**  paths.  However,  the  computation  is  simplified  by  the 

I identification  of  shipping  lanes  or  "isotropic  noise 

fields.”  In  either  case,  the  basis  is  a set  of  random 
ship  paths.  Each  is  u straight  line  (or,  e.g.  , a great 


L 


t 
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circle)  at  a random  range  and  makes  a random  angle 
with  respect  to  the  array  steering  angle. 


Ship  Path 


Steering  Direction 


‘ y 


1 

Array 


The  pairs  (y^,  tj/^)  are  independent  and  have  input  joint 

density  function  A single  ship  is  assumed  to 

travel  with  speed  v.  and  radiated  Intensity  SL.,  both 

J 

independent  of  time,  but  drawn  from  random  populations 
fo.v  different  ship  classes.  (Vj,  SLj)  are  independent 

pairs  with  density  function  fy.SL*  number 

of  ships,  K,  crossing  the  steering  direction  in  a unit 
time  interval  Is  assumed  to  be  a Poisson  variable  with 
mean  A (sbips/hour) ; 


PjK-kl  - e 


-X  y 
k! 


0,l,2,k • . 


(B-l> 
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This  is  the  key  to  the  shot-noise  process.  The  author 
of  the  model  describes  an  alternative  for  defining 
shipping  paths  in  terms  of  bearing  and  range  to  CPA. 

The  model  computations  are  simplified  as  the 
shipping  paths  become  more  structured.  A shipping  lane 
is  made  up  of  parallel  ship  paths  - so  that 

fy  * fy(y)  • (B-2) 

or,  alternatively,  a lane  may  be  defined  according  to 
constant  CPA  bearings,  with  random  CPA  ranges.  The  author 
also  defines  an  "isotropic  noise  field"  as  consisting 
of  many  paths  with  CPA  bearing  uniform  in  (0,  2it]  , and 
random  CPA  range.  Finally  a ship  scenario  might  consist 
of  several  Independent  lanes  plus,  perhaps,  an  indepen- 
dent isotropic  field. 

In  summary,  the  shipping  field  must  be  com- 
posed of  "straight"  paths  on  which  ships  arrive  according 
to  the  Poisson-rule  equation  (B-1).  The  grouping  of  paths 
into  lanes  or  isotropic  fields  makes  the  computations 
easier.  A ship's  speed  and  course  and  source  levels  do 
not  change  with  time,  but  are  chosen  randomly  for  each 
path  from  prespecified  probability  distributions.  The 
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source  intensity  does  not  depend  on  either  azimuthal  or 
vertical  angle.  All  of  the  ship  field  input  roust  be 
constructed  manually. 

2.2.4  Model  for  TL 

Transmission  loss  is  an  input  to  the  BTL  model, 
and  is  prescribed  as  a function  of  range  and  azimuth. 

It  is  treated  as  deterministic  and  there  is  no  roechanisn 
to  produce  random  fluctuations  of  the  TL.  The  more 
the  detail  of  the  TL,  the  more  complicated  is  the  calcu- 
lation of  the  statistics  of  noise.  The  model  treats 
neither  source  coupling  nor  receiver  response  in  terras 
of  multipath  or  multimode  transmission.  In  order  to 
derive  closed-form  results  for  special  cases,  Goldman 
has  assumed  an  A 4-  BlogR  rule  for  TL. 

2.2.5  * Receiver  Model 

Energy  from  each  noise  source  is  assumed  to 
arrive  from  the  source  direction  in  the  horizontal  plane 
as  a perfectly  coherent  plane  wave.  Hence,  the  array 
beamformer  is  modeled  with  the  usual  array-response 

■5  ^ 

Again,  see  Reference  B-3. 
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function  (beam  pattern);  it  is  a completely  determinis- 
tic input.  The  array  itself  is  assumed  to  be  located 
at  a fixed  geographic  position  and  depth. 

2.2.6  Details  of  the  Calculations 

For  each  ship  contributor,  the  noise  power  at 
the  array  is  given  in  general  by 


N^(t)  » h(t-tj^;  v^,  SL^,  y^,  (B-3) 


where  h is  simply 

SL^(t)-T^(t)-AG^(t), 

tj^  is  the  random  time  of  arrival  of  the  i-th  ship  at 
the  beam  axis.  Also, 


T^(t)  - T(R^(t),  0^(t)) 


AG^<t)  - AG(e^(t)), 
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where  R^(t)  Is  the  range  to  the  source,  and  0^(t)  is  the 
source  bearing  (relative  to  array  broadside),  both  func- 
tions of  Ivj^(t-t^),  y^l . 

Then  the  total  noise  intensity  at  time  t is  of 

form 


£h(t-t. ; 
i ^ 


v^,  SL^,  y^, 


(B-4) 


This  is  a "generalized  shot-noise  process"  (Reference  a-3), 
completely  characterized  by  the  ship  arrival  density  (X^), 
the  function  h,  and  the  distribution  functions  for  y^, 

Vj^,  and  SL^. 

i* 

Goldman  shows  that  the  ship  scenario  can  be 
equivalently  represented  as  a large  collection  of  ran- 
dom ship  paths  (with  a single  X)  or  in  terms  of  several 
shipping  lanes  plus  an  isotropic  field.  He  then  derives 
formulas  for  the  moments,  correlation  functions,  and 
characteristic  functions  for  N(t)  of  arbitrary  dimension 
(in  time).  For  example,  the  m-dimenslonal  characteristic 
function  is  given  by 


■•i--  r-»>'K.'C»»»W«'Sirp*»S^'?^^^SifW9teS**iS< 
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N(tj),N(t2),....N(tjjj) 


(B-5) 


where  E(  • ) denotes  expected  value  a;:id  j « 


Equations  of  similar  complexity  are  obtained  for  the  cumu- 
lants,  moments,  etc.  However,  for  special  cases,  simpli- 
fied formulas  are  obtained,  for  example,  if  there  is  one 
ship  path  only,  with  constant  values  of  v.  SL,  y,  li*.  then  for 
h(t)  - h(t-t^;  v,  SL,  y,  4»). 


E(N)  “ X 


f 

— iO 


h(t)dt. 


Covariance  (N)(t>  • X 


OD 


h(t)h(t+T)dt 


I 

1 

a ^ 


The  author  devotes  considerable  attention  to 
the  case  in  which  the  skewness 


C . SLI-Juf 
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is  small.  He  shows  in  Reference  (B~2)  that: 

for  a sequence  of  shot-noise  processes 
if  Cg  — ►O  and  the  correlation 
converges  to  p(t)  as  i->*'»,  then  the  Joint 
density  function  for  converges  to  a joint 
log-normal  density,  completely  character- 
ized by  the  mean,  variance  and  correlation 
function. 


The  validity  of  such  a proposition  yields  a very  simple 

model  for  noise  fluctuations.  The  condition  C —►0 

s 

is  shown  to  hold  for  a lane  whenever  one  of  the  following 
occurs : 


• X- 

• v--^o 

• R — ►<» 

Goldman  also  develops  (Ref.  B-2)  general  formulas 
for  the  Joint  characteristic  function  for  two  beams  of  one 
array  and  for  two  beams  ox  two  arrays. 

2.2,7  Output  and  Analysis 

The  model  output  consists  of  probability  den- 
sity functions  (of  arbitrary  order),  with  corresponding 
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moments  and  correlation  functions.  For  special  cases 
these  have  been  derived  in  closed  form;  otherwise  numeri- 
cal integration  or  Fourier  inversion  of  characteristic 
functions  is  required.  The  model  is  not  set  up  to  pro- 
duce sample  paths  (time  series  samples)  of  the  noise  process, 

2.2,8  Computer  Implementation 

Computer  routines  are  used  at  BTL  to  perform  the 
Fourier  inversions  (FFT),  numerical  integration,  etc.  The 
"model”  is  not  a general  production  computer  algorithm, 
so  that  each  case  and  type  of  output  must  be  considered. 

The  time  and  core  requirements  depend  on  the 
statistics  and  accuracy  desired.  For  the  UNI VAC  1108 
computer,  the  following  are  typical: 


- calculation  of  moments:  lOK  core  and  5 
seconds,  per  lane 


- calculation  of  one-dimensional  density  function: 
60K  core  and  5 minutes,  per  lane 


r 

I 


- calculation  of  two-dimensional  density  function: 
90K  core  and  12  minutes,  per  lane 


I 

I 
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2.2.9  Evaluation 

Output  of  the  BTL  model  has  been  compared  with 
surveillance-array  noise  data  and  the  agreement  has  been 
termed  ‘'excellent"  in  Reference  B-2.  In  particular,  the 
predicted  one-dimensional  density  functions  were  compared 
with  those  measured  over  a nine-day  period  for  two  steer- 
ing angles.  The  distribution  functions  were  within  0.5  dB 
over  the  range  0.02  to  0.97,  after  a shift  in  median  of 
3 dB. 


The  ship  arrival  assumption  was  also  tested 
against  extrapolated  observations  of  traffic.  The 
Poisson  model  was  accepted  at  a high  significance  level. 

2.2.10  Significant  Advantages  and  Disadvantages 

The  principal  contribution  that  the  BTL  work 
can  make  to  area  assessment  or  system  performance  analy- 
sis seems  to  be  the  analytic  formulas  derived  for  special 
"limiting"  cases.  These  can  be  used  not  only  to  predict 
the  important  beam  noise  fluctuation  properties,  but 
also  to  determine  the  environmental  and  shipping  para- 
meters which  influence  them.  As  examples  of  the  many 
results  derived  in  Reference  B-2,  consider 


2-30 


B - BTL 
Modwl 

• Noise  power  autocorrelation  depends  on 
“average  time  for  a ship  to  cross  the 
beam"  but  not  on  ship  density  (X),  width 
of  the  shipping  lane,  or  orientation  of 
paths. 

• As  the  skewness  of  the  one-<diroensional 
distribution  tends  to  zero,  the  noise 
process  approaches  a log-normal  process, 
depending  only  on  y,o  and  the  correlation 
function.  The  author  shows  that  the  con- 
dition is  satisfied  when  there  are  many 
ships  or  slow  ships  or  distant  ships.  The 
result  is  extended  to  the  multi-beam, 
multi-array  case. 

• Although  the  overall  ship  scenario  is 
very  important,  treatment  of  ship  speeds 
or  ranges  as  random  variables  (instead 
of  as  deterministic)  is  not  usually  ne- 
cessary . 


Even  though  these  results  do  not  take  into  account  TL 
fluctuation  or  the  details  of  the  array  response,  they 
are  indicative  of  the  power  of  the  approach. 


Of  all  the  Analytic  models  treated,  the  BTL 
model  is  the  most  general  in  the  sense  that  it  provides 
formulas  for  density  functions  of  all  diiwdnslons  and 
hence  a complete  description  of  the  level-crossing  pro- 
perties and  other  higher  order  statistics  of  the  noise. 
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It  can  also  treat  nearly  any  shipping  scenario  > subject 
only  to  the  Poisson  arrival  assumption.  On  the  latter 
point,  Goldman  has  examined  measured  ship-arrival  data 
and  found  the  hypothesis  to  be  acceptable  at  a high  sig- 
nificance level. 


Among  the  important  shortcomings  of  the  BTL 

model : 

• Except  for  simplified  cases,  the  calcula- 
tion of  characteristic  functions,  distri- 
bution functions,  correlation  functions, 
etc.  is  time-consuming  - especially  for 
the  higher  order  statistics. 

• The  preparation  of  shipping  lane  inputs  is 
difficult.  This  is  of  course  a problem  for 
any  of  the  models  which  use  more  than  ship 
densities. 

• The  asymptotic  limit  of  the  noise  process 
as  log-normal  is  an  extremely  powerful 
result.  However,  it  is  not  clear  that  it 
is  valid  in  interesting  cases  and  requires 
further  study. 

• The  key  assumption  about  Poisson  ship 
arrivals  needs  further  study . both  in  the 
sensitivity  of  the  analytic  results  to  it 
and  in  its  validity  as  a property  of  shipping 
traffic  behavior. 
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In  addition  to  these  is  the  ensembling  problem  discussed 
in  Section  1.  The  density  functions  produced  by  the  model 
are  derived  from  the  ensemble  over  all  ship  locations. 

Since  the  shot  process  is  (in  most  cases)  ergodic,  the 
ensemble  statistics  are  the  same  as  time-average  statistics, 
over  time  periods  very  long  compared  to  the  decorrelation 
ti*^.  Hence,  the  density  function  for  a 12  or  24  hour 
time  period  will  usually  not  be  represented  by  the  BTL 
model  prediction  of  the  one-dimensional  density.  This  is 
not  to  say  that  the  model  is  incapable  of  predicting  the 
properties  of  short  sample  paths;  the  BTL  formulation  can 
in  theory  provide  such  data,  but  the  calculations  are 
expected  to  be  very  tedious. 
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2.3  BBN  MODEL 

2.3.1  Background 

Name:  (BBN  Noise  Model)  (BBN) 

Developer:  Bolt.  Beranek  and  Newman,  Inc., 

John  I.  Mahler,  Francis  J.  M.  Sullivan, 
Magnus  Moll 

Sponsor:  Office  of  Naval  Research,  Code  431; 

Naval  Electronics  Systems  Command. 

Code  320 

Previous  Application: 

Published  Documentation:  References  C-1,  C-2. 

2.3.2  General  Approach 

The  BBN  oiodel  uses  an  Analytic  approach,  based 
on  Equation  (1-1).  It  involves  the  direct  calculation  of 
one-  and  two-dimensional  (in  time)  probability  density 
functions  by  inverse  Fourier  transform  of  the  characteristic 
functions.  The  key  to  the  approach  is  the  construction  of 
the  characteristic  function  when  ship  traffic  and  source 
levels  are  modeled  as  random  variables.  (General  shipping 
scenarios  are  treated  and  a special  model  for  radiated 
noise  encompassing  the  discrete  and  con  inuous  parts  ox 
the  spectrum  is  used  (Ref.  C-3). 
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In  the  current  version  of  the  model  only  first- 
order  statistics  are  computed.  The  computer  implementation 
of  the  method  for  second-order  densities  has  not  been  com- 
pleted. Furthermore,  only  "main  beam"  noise,  i.e.,  noise 
from  a fixed  azimuthal  sector,  is  provided  at  present. 
Transmission  loss  and  array  response  are  deterministic 
inputs  to  the  routine. 

The  BBN  model  is  similar  to  the  BTL  model  in  that 
both  calculate  density  functions  directly  from  characteristic 
functions.  The  primary  difference  is  that  BTL  assumes  that 
the  ships  travel  according  to  a Poisson  rule  (which  in 

turn  simplifies  the  calculation  of  the  characteristic  func- 

♦ 

tion)  while  BBN  allows  a general  ship  scenario. 

The  time  for  the  computer  calculation  of  the  charac 
teristic  function  and  Fourier  inversion  for  the  one-dimen- 
sional density  (with  £0-60  dB  dynamic  range)  is  about  one 
CDC-6400  minute  ($10).  Two-dimensional  densities  are 
expected  to  utilize  much  more  computer  time. 


♦The  version  of  the  BBN  model  most  recently  implemented  is 
described  in  Reference  C-2.  It  is  assumed  there  that  the 
number  of  ships  in  the  "urea  of  interest"  Is  a Poisson 
variable,  as  in  the  UST  model, 
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2.3.3  Model  for  Ships 

Ships  are  confined  to  a bounded  domain,  the  basin, 
as  in  the  two  models  described  above.  Merchant  ship  traffic 
is  assumed  to  travel  within  "route  envelopes"  or  "lanes." 

For  the  basic  time  interval  of  interest,  say  0 s t < T, 
the  model  need  account  only  for  the  ships  which  pass  through 
the  azimuthal  sector  under  consideration  (the  beam).  Hence, 
for  each  lane  the  "Area  of  Interest"  in  which  ships  are  de- 
fined looks  like: 


Receiver 


so  that  it  contains  ships  which  could  reach  the  beam  with- 
in 0 i t < T. 


i 
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A shipping  scenario  is  made  up  of  (a)  individual 
ships  of  (b)  several  types  or  classes,  traveling  on  (c) 
tracks  within  (d)  routes  or  lanes. 


Routes : The  scenario  consists  of  M routes. 

Tracks : A route  is  made  up  of  straight-line 
tracks.  Track  bearings  are  inputs. 

All  tracks  within  a route  are  equispaced 
in  range  when  viewed  along  the  center  of 
the  array’s  main  beam.  The  number  of 
tracks  within  a route  (track  density) 
is  chosen  so  that  the  resulting  resolu- 
tion is  appropriate  to  the  noise  pro- 
blem under  study. 


Types : Each  ship  is  assumed  to  belong  to  one 
of  N types  or  classes.  Then,  for  the 
i-th  route  there  is  a random  variable, 

Jf  j ’ determines  the  number  of  ships 

of  the  J-th  type  in  the  area  of  interest.* 
There  is  also  a constant,  deterministic 
speed,  Vj,  and  a source-intensity  random 
variable,  SL^  (based  oh  Ref.  C-3):  these 
are  the  same  for  all  ships  of  type  j. 


The  version  of  the  model  described  in  Reference  C-2  assumes 
thac.  Jh4  is  a Poisson  variable.  The  assumption  is  based 
on  a simple  model  of  traffic  which  leads  to  a binomial 
distribution,  which  in  turn  is  approximated  by  the  Poisson 
distribution  for  small  segments  of  the  route  associated 
with  small  azimuth  sectors. 
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Individual  Initial  ship  positions  within  the 
Ships:  region  of  interest  are  determined 
from  the  random  variables: 


“ coordinate  in  the  direction  of 
travel  for  the  k-th  ship  of 
type  J on  route  1 


“ coordinate  normal  to  the  direc- 

^ijk 

tion  of  travel  for  the  k-th 
ship,  etc. 


[ 

I 

I 
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These  variables  determine  positions  on 
a track  and  across  the  tracks  for  a route. 
They  are  independent  from  ship  to  ship, 
but  are  identically  distributed  as  a 
function  of  k.  The  density  functions 
for  G and  Q are  inputs  to  the  model . The 
number  of  ships  on  a route  and  of  a type 
is  mentioned  above. 

Notice  that  Q is  the  location  at  time  of 
CPA  and  G is  a multiple  of  the  ship  velo- 
city. Since  ships  travel  at  constant 
speed  (Vj  for  type  J),  the  positional 
coordinates  at  time  t are 

Q(t)  - Q. 

G(t)  - G + Vj.t 


li 
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2.3.4  Model  for  TL 

Transmission  loss  is  an  input  to  the  BBN  model. 
It  is  a deterministic  function  of  range  appropriate  to 
the  azimuthal  sector  under  consideration.  The  source  and 
receiver  models  do  not  depend  on  vertical  transmission 
or  vertical  arrival  angle. 

2.3.5  Receiver  Model 

Energy  from  each  noise  source  is  assumed  to 
arrive  from  the  source  direction  in  the  horizontal  plane 
as  a perfectly  coherent  plane  wave.  The  "array  response" 
is  for  an  ideal  beam:  unit  response  over  the  prescribed 
azimuthal  sector  and  zero  response  elsewhere.  Energy 
from  different  sources  is  added  on  a random  phase  basis. 

2.3.6  Details  of  the  Calculation 

Begin  with  equation  (1-1): 

j(^) 

N(t)  - SLj(t)’Tj(t)-AGj(t).  (1-1) 

For  the  ships  in  the  sector  of  interest,  AG  h 1.  The 
sura  can  be  written  in  terms  of  routes,  types,  and  indivi- 
dual ships  after  we  define: 
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• * source  intensity  of  the  k-th  ship 

of  the  j-th  type  on  route  i , a random 
variable,  constant  in  time. 


^^ijk’  ^ijk^ 


the  coordinates  at  time 
t«0  of  the  k-th  ship,  random  variables 
constant  in  time  as  described  in  2.3.3. 


• T(G,  Q)  » the  transmission  ratio  (inten- 
sity loss)  for  a source  at  coordinates 
(G,  Q).  It  is  a deterministic,  time- 
independent  variable. 

• M ■ number  of  routes 


• N « number  of  ship  types 


• '^ij  * *'^***^®**  ships  of  type  j in  the  area 
of  interest  for  the  i-th  route. 


Then  (1-1)  at  time  t«0  becomes 

M N “^ij 

N(0)  - £ E S SL 
i"l  J»1  k»l 


Since  ships  of  the  same  type  travel  at  constant  speed, 

Vj  the  noise  intensity  at  time  t is  given  by 

M N '^ij 

■ S S k?i  «ljk)  ‘C-2) 
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These  two  equations  (C-1,  C-2)  are  the  basis  for 
the  BBN  approach.  Notice  that  the  random  variable  N(t) 
is  a function  of  the  random  variables  SL,  G,  Q,  and  J, 
and  of  the  deterministic  variables  v,  t,  M and  N.  Hence, 
the  density  functions  and  statistics  of  N are  derived 
from  ensembling  of  SL,  G,  Q,  and  J. 


Ref.  (C-1)  mentions  that  the  calculation  of 
such  densities  can  be  performed  directly  with  convolutions. 
This  might  proceed  as  follows  for  the  one-dimensional  den- 
sity : 

• For  each  i,  j and  t fixed,  calculate 

«ijk>  1 <=} 

P(SLijk  <_  D) 

P(Jij  1 E)  (C-3) 


These  are  independent  of  k. 


• For  a single  ship  of  type  J on  the  i-th 
route,  find 


P 


which  is  again  independent  of  k.  If  SL 
and  T are  in  dB's,  then  this  is  a convolu- 
tion. 
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• Given  (C-3)  and  (C-4),  for  a single  (i.j), 
compute 


which  again  involves  multiple  convolutions 
for  the  first  term . 


• Finally,  N(t)  is  simply  a double  sum  of 
the  variables  given  in  (C-5)  - so  that 
M*N  additional  convolutions  yield  the 
answer. 

Calculation  of  the  joint  density  for  (N(tj^),  N(t2))  is 
similar,  but  involves  the  joint  distribution  of  form 

p|T(G+vtj.  Q)  < C,  TCG+vtg,  Q)  < d|. 


For  implementation  on  a digital  computer,  the 
use  of  characteristic  functions  and  Fourier  transforms 
is  often  a better  approach  than  that  described  above. 

The  present  BBN  model  does  indeed  employ  such  a method. 
Reference  (C-1)  derives  the  formula  for  the  characterise 
tic  function  corresponding  to  the  joint  density  of 
(N(0).  N<t)); 
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♦(a, 6) 


M N 

n n 


gij(w) 


•t  (w)dw 

-I 


(C-6) 


where 


I 


and 


is  the  characteristic  function  for  SL 


T 


i 

i. 


i. 


g. . is  the  density  function  (for  all  k)  of  the 

^ w 

variable 


aT(G 


ijk’ 


^ijk^ 


+ 6T(G 


ijk 


^ijk^’ 


The  latter  density  is  found  from  its  characteristic 
functions. 


Since  <Ka,6)  ® Ejexp(iaN(0)  + iBN(t))j,  the 
characteristic  function  for  N(t)  is  simply  ♦(0,6).  At 
present,  the  BBN  program  computes  only  the  one-dimensional 
density  by  applying  an  FFT  to  the  characteristic  function 
of  N(t). 
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2.3.7  Output  and  Analysis 

The  model  output  consists  of  the  one-dimensional, 
discretized  probability  density  function  for  noise  inten- 
sity. The  dynamic  range  can  be  50-60  dB.  All  first  order 
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moments  can  then  be  found  directly.  The  model  does  not 
produce  sample  paths  or  any  higher  order  statistics  at 
present;  however,  the  two-dimensional  (in  time)  density 
and  corresponding  second  order  statistics  (e.g.,  auto- 
correlation function)  will  be  the  outputs  of  a future  ver- 
sion. 

2.3.8  Computer  Implementation 

The  BBN  model  is  coded  in  FORTRAN  and  imple- 
mented on  the  CDC-6000-series  computer.  For  a typical 
run,  to  determine  the  one-dimensional  density,  the 
calculation  requires  a 10-second  FFT  and  a total  time  of 
about  one  minute  on  the  CDC  6400  ($10).  The  amount  of 
time  for  computing  the  two-dimensional  density  is  expected 
to  be  large. 

2.3.9  Evaluation 

Comparisons  of  model  output  with  experimental 
data  (from  LAMBDA)  has  been  proposed;  but  has  not  yet 
begun . 

2.3.10  Significant  Advantages  and  Disadvantages 

Like  the  BTL  model,  the  BBN  approach  provides 
an  analytic  determination  of  the  noise  density  function. 
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Carefully  chosen  cases  can  be  studied  to  determine  the 
bounds  and  driving  parameters  for  the  noise  fluctuations. 

The  use  of  general  source-intensity  and  ship 
traffic  models  are  special  features  of  the  BBN  program. 

The  model  requires  the  manual  construction  of 
ship  scenarios  and  input  of  TL  and  other  data.  The  most 
significant  shortcomings  are  the  lack  of  higher  order 
statistics  and  the  ensembling  problem  discussed  in 
Section  1 and  in  connection  with  the  USX  and  BTL  models. 
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2.4  WAGNER  MODEL 

2.4.1  Background 

Name:  (Wagner  Associates'  Noise  Model)  (WAGNER) 

Developer:  Daniel  H.  Wagner,  Associates, 

Bernard  J . McCabe 

Sponsor:  Office  of  Naval  Research,  Code  431 
Previous  Applications:  Sonobuoy  Problems 
Published  Documentation:  Reference  D-1 

2.4.2  General  Approach 

The  WAGNER  model  was  developed  for  the  sonobuoy 
problem  and  hence  does  not  explicitly  predict  beam-noise 
statistics.  It  is  included  in  this  survey  because  it 
provides  analytical  results  which  supplement  those  of 
the  other  models  and  because  it  could  potentially  be  modi- 
fied to  produce  array  noise  data.  Furthermore,  the  model 
considers  not  only  the  temporal  correlation  of  ship-generated 
noise,  but  also  the  spatial  correlation  of  average  intensi- 
ties for  distributed  sensors. 

The  WAGNER  model  actually  includes  three  differ- 
ent approaches:  two  of  them  are  Analytic  and  complementary. 
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while  the  third  is  Brute  Force  and  based  on  the  analytic 
formulation.  Unlike  the  previous  models,  the  analytic 
approaches  are  not  computer-coded  to  yield  density  func- 
tion, etc.,  but  rather  are  used  to  provide  characteristic 
functions  and  approximate  results  which  in  turn  can  be 
used  to  obtain  simplified  models  and  insight  into  the 
noise  process.  The  Brute-Force  model  will  be  described 
here  only  in  terms  of  its  analytic  basis. 

The  WAGNER  approach  is  based  on  equation  (1-1), 
but  is  distinguished  in  the  treatment  of  the  shipping 
scenario  (as  are  the  other  models  reviewed  above). 

2.4.3  Model  for  Ships 

The  two  WAGNER  ship  models  are  termed:  ’’Poisson” 
and  "Bounded”  noise  processes. 

Poisson:  At  t»0  the  noise  sources  are  distri- 
buted according  to  a 2-dimenslonal 
(spatial)  Poisson  process  with  inten- 
sity (density): 

X • mean  number  of  ships  per  unit  area. 

This  means  that  the  number  J of  ships 
in  a region  of  area  A is  a Poisson 
variable  with  parameter  XA: 

P<J-k)  - e'*^*  ' (XA)‘‘/k:.  (D-i) 
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Moreover,  the  numbers  of  ships  in  two 
disjoint  regions  are  independent  Poisson 
variables;  and  finally  the  ships  in  a 
region  are  uniformly  and  independently 
distributed  over  that  region. 

Bounded:  At  time  t=»0,  a fixed  finite  number  of 
ships  are  selected  independently  from 
a uniform  distribution  over  a bounded 
region. 

Notice  that  the  Poisson  process  is  similar  to  the  US I 
model  in  that  the  ships  are  uniformly  distributed  in  a 
region  aind  that  the  ship  count  is  Poisson.  On  the  other 
hand,  it  assumes  an  unbounded  domain  and  constant  density 
over  that  domain.  The  author  of  the  WAGNER  model  notes 
in  Reference  D-1  that  the  Poisson  process  is  useful  for 
obtaining  theoretical  results  while  the  Bounded  process 
is  needed  to  simulate  real  environments. 

For  either  process,  initial  ship  courses  are 
uniform  random  variables  on  (0,360),  san^led  independently 
for  each  ship.  Likewise,  there  is  an  input  ship-speed 
distribution  and  a source- intensity  distribution  from 
which  satqples  are  taken  independently  for  each  ship.  Once 
chosen,  those  variables  remain  ^..onstant  in  time,  with  the 
exception  that  in  the  Bounded  process  ships  reflect  from 
the  boundary,  with  angle  of  incidence  equal  to  angle 
of  reflection.  All  ship  sailing  is  "rectangular.'  i.e.. 
as  if  on  a flat  earth. 
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2.4.4  Model  for  Transmission  Loss 

Transmisaion  loss  is  a deterministic  input,  de- 
pending only  on  range.  No  aaimuthal  dependence  is  allowed, 
nor  are  the  effects  of  vortical  arrival  structure,  co- 
herence, etc.  included. 

2.4.5  Rece Ivor  Modei 

The  receiver  is  assumed  to  be  an  omnidirectional 
hydrophone,  fixed  in  location.  Its  response  is  a deter- 
ministic input,  which,  together  with  the  TL  and  sourct' 
level,  accounts  for  temporal  processing.  For  the  WAGNER 
analytic  approach,  the  ensembling  problem  discussed  above 
is  prOvsent. 


The  results  for  spatial  correlation  are  for  two 
separated  oimit  phones.  The  correlation  is  in  intensity, 

so  that  relative  phases  are  luit  required. 

2.4. G Details  of  ttu^  Calculation 

The  calculalton  t)f  the  noi.se  proct'ss  Is  based 
on  equation  (l-l): 

N(t)  ="  £ SI.  -T  (t), 

.K.I  '' 

wheri'  .)  is  a countable  index  set.  Refc'rence  D-l  first 
shows  that  for  t fixed  and  T(rj)  equal  to  the  t ransmissii>n 
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ratio  from  the  receiver  to  a 

Poisson  process  at  range  r ^ 

J ‘ 


for  N(t)  - 


J-1 


sample  from  the  spatial 
the  characteristic  function 


00  00 

$(w)  » exp  1 271  xy* y*  (e^^*''^^*'^-l)rdrdF(z)| , (D-2) 


o 0 

where  F is  the  conunon  distribution  function  for  SL 


‘j 


then 


It  follows  from  (D-2)  that  when 


00 

./rT 


T » / rT(r)  dr  < », 


and 


E(N)  - 27iXE(SL)*T  < « 

00 

= 2tiXe(SL^)  y rT^(r)dr  < 


(D-3) 


(D-4) 

(D-5) 


The  formulas  (D-2),  (D-4),  and  (D-5)  are  quite  useful  in 
estimating  first-order  statistics  for  N(t). 

McCabe  (Ref.  D-1)  also  shows  that  N(t)  is  sta- 
tionary for  the  Bounded  process  on  a rectangle  and  for 
the  Poisson  process.  He  then  proposes  several  theses: 

(a)  Although  no  analytic  method  has  yielded  an 
autocorrelation  function  for  N(t),  Monte- 
Carlo  samples  for  a smooth  TL  (with  conver- 
gence zone  structure)  suggest  that  noise 
levels  N(t)  (in  dB)  have  a nearly  exponential 
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autocovariance  function,  of  form 

C(t)  ■ o*  exp(-t/T).  (D-6) 

(b)  In  the  san^ile  replications,  the  "decorrela- 
tion" times,  T,  for  (D-6)  agree  well  with 

1/T  = (1.5)  V (D-7) 

where  v is  average  ship  speed  and  X is  the 
Poisson  ship  density  (ships  per  unit  area). 
This  formula  can  be  explained  by  the  argu- 
ment that  the  decorrelation  time  should  be 
related  to  the  waiting  time  for  the  "nearest- 
neighbor"  ship  to  be  replaced.  Nearest 
neighbor  changes  occur  every  1/vn/x  time 
units  in  the  case  that  ships  are  uniformly 
distributed  and  move  at  average  speed  v. 

(c)  The  noise  level  N(t)  (in  dB’s)  is  approxi- 
mately Gaussian  (all  dimensions).  The  author 
quotes  a central-limit  theorem  of  Marlow 
(Ref.  D-2). 

(d)  Finally,  since  N(t)  is  stationary,  nearly 
Gaussian,  and  has  an  autocovariance  func- 
tion which  is  approximately  exponential, 
the  author  concludes  that  the  process  is 
approximately  Gauss-Markov  (by  Doob’s 
Theorem,  Ref.  D-3). 

If  the  last  proposition  were  valid,  then  the  noise  time 
series  could  be  simulated  very  easily  and  inexpensively 
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(compared  to  the  Brute  Force  methods  or  other  Analytic 
methods).  Like  the  BTL  asymptotic  limit,  this  approach 
merits  further  investigation,  especially  for  possible 
extension  to  the  beam-noise  problem. 

Reference  D-1  also  deals  with  the  spatial  corre- 
lation of  ambient  noise.  A formula  for  the  (0  time  lag) 
cross-correlation  of  N(t)  at  two  separated  points  is  de- 
rived for  the  Poisson  process.  It  does  not  depend  on  the 
values  of  X (density  of  ships)  or  the  source-level  distri- 
bution, but  is  driven  by  the  TL.  For  the  noise  levels 
(N(t)),  it  is  shown  that  under  certain  conditions  the 
same  conclusion  holds  (correlation  depends  only  on  TL). 
Monte-Carlo  simulations  suggest  that  significant  variations 
in  TL  (e.g.,  convergence  zones)  cause  the  spatial  correla- 
tion to  tend  to  zero  quickly  (i.e.,  within  a short  distance). 

The  WAGNER  approach  includes  a method  for  simu- 
lating the  noise  for  a field  of  omni  sensors.  It  involves 
the  construction  of  a set  of  Gauss-Markov  processes  which 
have  the  proper  spatial  correlation. 

2.4.7  Output  and  Analysis 

The  general  approach  provides  formulas  for  the 
one-dimensional  characteristic  function,  the  first  two 
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moments,  and  the  autocorrelation  function  for  the  omni 
noise  intensity.  If  the  Gauss-Markov  approximation  is 
valid,  then  formulas  for  higher-order  statistics  and 
efficient  simulation  routines  for  the  noise  levels  are 
available. 

2.4.8  Computer  Implementation 

These  are  routines  simulating  the  Gauss-Markov 
process  and  the  Brute-Force  method  (not  covered  here). 
Moreover,  a computer  program  for  calculating  the  spatial 
correlation  is  available.  A typical  run  with  1000  Poisson 
samples  and  correlation  calculated  at  5 nm  increments 
from  0 to  100  nm  costs  about  $50. 

2.4.9  Evaluation 

The  results  have  not  been  compared  with  measure- 
ment data. 


2.4.10  Significant  Advantages  and  Disadvantages 

As  discussed  in  2.4.1,  the  WAGNER  model (s)  has 
limited  direct  applicability  to  the  beam-noise  problem. 

It  assumes  an  omni  sensor  and  a homogeneous  ship  distri- 
bution. There  are  however  useful  analytic  results  which, 
as  in  the  case  of  the  BTL  model,  should  not  be  ignored 
in  the  synthesis  of  an  approach  to  predicting  beam-noise 


statistics. 
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2.5  NABTAM 

2.5.1  Background 

Name:  (Narrow  Beam  Tewed  Array  Model)  (NABTAM) 

Developer  : Raff  Associates:  W.  Galati, 

£.  Moses,  R.  Jennette.  (Operations 
Research,  Inc.,  Underwater  Systems, 
Inc. ) 

Sponsor  : Office  of  Naval  Research  (Code  431), 
LRAPP,  NORDA 

Previous  Applications:  Studies  for  LRAPP, 

DDR&E,  and  ONR  431 

Documentation:  In  Preparation  (Reference  E-1) 

2.5.2  General  Approach 

NABTAM  Is  a Brute-Force  model  designed  to  pre- 
dict the  response  of  a horizontal  line  array  to  wind- 
generated  noise,  surface  ships,  and  a target  at  a single 
frequency  in  the  range  of  10-1000  Hz.  The  wlnd-nolse  and 
target  signal  do  not  change  In  time.  The  time-dependent 
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surface-ship  noise  is  calculated  from  equation  (1-1),  and 
the  remainder  of  this  description  of  NABTAM  concentrates 
on  that  aspect  of  the  model. 

Some  attributes  which  distinguish  NABTAM  from 
most  of  the  other  Brute-Force  models  are: 

(a)  The  transmission  loss  from  source  to  re- 
ceiver is  calculated  within  the  program. 

It  is  a ray-trace  routine  for  a range- 
independent  environment . 

(b)  A version  of  the  program  (maintained  by 
ORI)  performs  "near-field"  corrections, 
i.e.,  it  calculates  array  response  for 
cylindrical  wavefronts  from  nearby  sources. 

(c)  The  TL  and  receiver  models  take  into  account 
the  vertical  arrival  structure  of  the  noise 
intensity. 

(d)  Beam  noise  for  multiple  array  locations, 
depths,  and  orientations  can  be  calcu- 
lated in  one  model  run  (i.e.,  for  the  same 
ship  field).  Beam  patterns  can  be  calcu- 
lated internally. 


The  model  is  designed  to  generate  a single 
realization  of  the  noise  time  series  for  one  or  several 
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array  beams.  (Monte  Carlo  simulations  are  not  generally 
performed,  nor  is  code  structured  to  produce  the  ensemble 
statistics).  Initial  ship  locations  are  inputs.  Ship 
tracks  arc  chosen  from  probability  distributions,  but 
remain  constant  for  the  duration  of  the  replication.  All 
other  parameters  (TL,  receiver  response,  source  levels, 
array  location)  are  constant  in  time  and  deterministic. 

MABTAM  is  programmed  in  FORTRAN  IV  and  Installed 
on  a number  of  computers,  including  the  CDC  6600  at  Eglin 
AFB. 


2.5.3  Model  for  Ships 

All  surface  ships  are  assumed  to  have  the  same 
source  level , as  determined  from  the  Ross  and  Alvarez 
"normal  merchant  ship"  (Ref.  E-2).  This  level  is  con- 
stant in  time,  and  depends  only  on  acoustic  frequency. 
The  source  intensity  does  not  depend  on  ship  speed  or 
on  transmission  angle. 

Initial  ship  locations  are  a deterministic 
input.  There  is  no  real  limit  to  the  number  of  ships 
allowed,  but  no  new  ships  can  be  created  after  the 
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initial  time  step.  The  initial  ship  courses  and  speeds 
can  be  user-specified  or  they  can  be  drawn  from  random 
populations.  In  the  latter  case,  the  parameters  are 
determined  independently  for  each  ship,  with  the  course 
taken  from  a uniform  distribution  over  a specified  set 
of  angles  and  the  speed  approximately  normal  (mean  and 
variance  are  inputs).  Once  established,  ships  travel  on 
tracks  determined  by  d(Latitude)/dt  and  d(Longltude)/dt 
constants.  Hence,  even  speed  is  not  necessarily  constant. 
If  a ship  leaves  the  basin  area,  it  is  not  replaced. 

There  is  no  mechanism  for  doing  replications, 
except  to  rerun  the  model. 

2.5.4  Model  for  TL 

NABTAM  is  the  only  noise  model  reviewed  here 
which  has  its  own  internal  TL  routine.  The  program  is 
structured  in  such  a way  that  multiple  receiver  depths 
can  be  treated  in  one  run  and  vertical  arrival  structure 
is  taken  into  account.  Hence,  the  TL  as  a function  of 
range,  depth,  and  vertical  arrival  angle  is  required. 

The  program  actually  precalculates  the  ray-trace  para- 
meters as  functions  of  vertical  angle  and  can  thus 
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generate  TL  for  each  source  location  at  each  time  step 
(but  for  a single  frequency). 

The  TL  model  is  based  on  geometric  acoustics, 
with  no  corrections  for  diffraction.  It  assumes  a range- 
independent  environment:  sound  speed,  water  depth,  and 
boundary  los.ses  do  not  change  in  range.  The  sound  speed 
profile  is  linearly  segmented  in  depth.  Rays  are  traced 
according  to  Snell's  Law  from  the  receiver  depth  in  1° 
vertical-angle  steps  for  one  cycle.  If  the  1°  increments 
are  not  sufficient  (determined  in  the  intensity  calculation), 
then  additional  rays  are  traced  in  the  angular  regions  when 
they  are  needed.  Rays  are  classified  according  to  ten 
families  (e.g. , RR,  RSR,...)  and  the  range  at  which  the 
source  depth  is  crossed  is  stored.  The  spreading  loss  is 
calculated  for  each  path  from  the  usual  expression 


-lOlog 


cos9 


» 


where  R is  range,  0 is  receiver  angle,  and  0^  is  source 
angle.  The  value  of  dR/dO  is  estimated  from  an  Interpola- 
tion between  rays  of  the  same  family  which  bracket  the 
source  range  at  the  source  depth.  This  is  the  approach 
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used  in  many  of  the  early  ray-trace  models  (e.g. , FAST 
NISSM  (SHARPS  II)  or  RP-70),  but  has  much  less  precision. 

The  intensity  becomes  infinite  at  caustics  (where  dR/d6  = 0), 
and  the  model  uses  preselected  truncation.  All  paths  are 
added  on  a random-phase  basis  so  that  detailed  multipath 
interference  and  surface- image  interference  are  not  pre- 
dicted. Individual  arrivals  at  the  receiver  are  collected 
and  summed  in  vertical-angle  bins  so  that  TL  as  a function 
of  vertical  angle  for  the  range  to  the  source  is  available 
for  the  calculation  of  array  response. 

The  TL  is  deterministic  and  time-invarient . 

2.5.5  Receiver  Model 

The  program  can  calculate  the  array-response  time 
series  for  multiple  receiver  locations  and  orientations, 
using  the  same  source  field,  in  one  program  execution. 
Receiver  locations,  depths  and  beam  patterns  are  determin- 
istic program  inputs  which  r^ain  constant  once  specified. 
Each  receiver  is  an  array  whose  response  to  plane-wave* 
arrivals  is  characterized  in  terms  of  beam  patterns, 
accounting  for  vertical  and  azimuthal  angles,  main  beams 
and  sidelobes.  All  energy  is  assumed  to  arrive  as  plane 

*The  ORI  version  of  NABTAM  models  the  near-field  response 
with  cylindrical  wavefronts. 
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waves,  distributed  in  vertical  angle,  but  perfectly  coherent 
and  undistorted.  The  arrivals  from  a given  source  (as 
calculated  from  the  TL  routine)  are  sorted  according  to 
vertical  angle  and  ’’convolved'’  with  the  beam  pattern. 

For  user  convenience,  NABTAM  has  routines  which 
calculate  beam  patterns  for  uniformly-spaced,  horizontal 
line  arrays. 

2.5.6  Details  of  the  Calculation 

Once  the  receiver  parameters,  shipping  field,  envi- 
ronmental inputs  and  array  response  functions  are  specified, 
the  noise  is  calculated  at  discrete  time  steps  according 
to  equation  (1-1).  Readings  are  taken  at  pre-spec ified  time 
points  for  a selectable  time  span.  The  ship  positions  are 
tracked  automatically  in  time  and  are  the  single  cause  for 
the  time  dependence  of  the  beamformer  output. 

The  flow  of  the  computer  program  is  shown  in  Figure 

E-1. 


Ve  note  that  wind-generated  noise  is  a time- independent 
contribution  added  to  the  ship-noise  field.  The  model  treats 
wind  sources  as  continuously  distributed,  uncorrelated  point 
sources  at  the  ocean  surface.  Provision  is  made  to  divide 
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the  surface  into  angular  sectors,  each  specified  radial 

extent.  The  contributions  in  each  horiz<.  sector  are 
computed  in  one-degree  vertical  stcpo,  using  a technique 
based  on  an  approach  suggested  by  Talham  (Ref.  E-3),  which 
permits  the  contribution  at  each  vertical  angle  to  be  com- 
puted from  a closed  form  expression.  The  resulting  noise 
intensity  is  thus  a function  of  both  horizontal  and  vertical 
angle. 


Finally,  the  N.ABTAM  model  is  structured  to  produce 
array  response  to  fixed  targets,  as  well  as  beam-noise  time 
series.  In  the  former  application,  the  usual  program  opera- 
tion is  to  calculate  the  array  response  for  a set  of  target 
ranges  and  bearings,  specified  relative  to  the  array.  For 
the  latter  case,  the  receiver  location  and  depth  can  be 
varied  but  the  noise  source  (ship)  locations  at  each  time 
step  are  fixed.  Combining  the  two  sets  of  results  permits 
estimates  of  signal  to  noise  ratio  as  a function  of  IxUh 
array  and  tat'get  pv>»ition.  while  avoiding  the  massive  number 
of  runs  required  to  handle  each  target-receiver  p^^sition 
combination  as  an  individual  model  calculation. 

2.5.7  Output,  and  Analysis 

The  output  of  the  NABTAM  models  consists  of  a beam- 
noise  time  series  for  each  beam  pattern  (e.g. . for  several 
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steering  directions),  for  a single  frequency,  as  a function 
receiver  location.  No  other  analysis  is  performed,  although 
a time-series  package  (e.g.  , as  in  DSBN)  could  clearly  be 
used  tc  generate  statistics  of  the  model  output.  As  noted 
above,  the  model  can  also  predict  wind-noise  and  target- 
signal  level?  (both  independent  of  time)  at  the  array  output, 
and  combine  them  with  the  ship  noise  to  yield  signal-to- 
noise  ratios. 

2.5.8  Computer  Implementation 

The  NaBTA’  Kiel  is  programmed  in  FORTRAN  IV.  It 
has  been  run  on  GE-635,  CDC  6000-series,  and  other  computers. 
The  shipping  parameters  (and  beam  patterns,  if  '.ot  for 
specific  line-array  types)  are  user  inputs  and  are  not  auto- 
matically generated.  The  program  is  overlaid  and  the  core 
storage  requirement  is  about  lOK  words.  If  it  were  not 
overlaid,  about  13K  words  would  be  needed. 

Costs  depend  on  the  number  of  ships,  beam  patterns, 
time  steps,  time  span,  etc.  As  an  example,  one  time  step  for 
a large  number  of  ships  and  five  beams  took  9 CDC-6400  .seconds 
($0.75).  A calculation  of  beam  noise  for  array  locations  and 
hundreds  of  surface  ships  for  a limited  number  of  time  steps 
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can  cost  $100  or  more.  If  the  program  were  not  overlaid, 
these  costs  would  be  reduced.  Also  note  that  NABTAM  calcu- 
lates TL  internally,  an  important  consideration  in  comparing 
noise-model  computer  costs. 

2.5.9  Evaluation 

There  has  been  no  foi'mal  evaluation  of  the  NABTAM 
model.  Comparisons  were  made,  but  not  documented,  several 
years  ago  with  Ionian  Basin  (Med)  and  LRAPP  omnidirectional 
noise  data.  The  levels  predicted  were  lower  than  measured 
by  several  dB,  but  trends  in  time  were  reportedly  simulated 
accurately. 

2.5.10  Sign! f lean t Advan la^es  and  Disadvantages 

Some  of  the  attribute's  which  distinguish  NABTAM 

from  other  Brute-Force  modt>ls  are; 


(a)  The  three-dimensional  distribution  of  energy 
w \ t h ang  I e a t t lu*  rt'ce  i ver  is  cal  cu  I a t ed  for 
each  m>ise  source.  Hence,  the  time-dependent 
noise  output  t.>f  a volumetric  or  vertical 
lint'  array  cun  bi'  predicted.  Most  impm'tunt 
for  this  revit'w,  tht'  bi'am  output  iM'  a hori- 
v.ontal  line  array,  whiut  steered  to  off- 
broadside  directions,  can  be  computed  with 
.multipath  "beam-splitting"  properly  accounted 
for . 
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(b)  The  program  is  structured  to  treat  multiple 
receiver  locations  and  response  patterns 
with  little  computational  redundancy. 

(c)  Array  output  for  wind  noise  and  targets  can 
be  calculated  in  the  same  execution  as  that 
for  ship  noise. 


Among  the  limitations  intrinsic  to  NABTAM  are: 


(a)  The  surface  ship  model  will  eventually  break 
down  as  ships  leave  the  basin  of  interest; 
the  ship  paths  and  speeds  are  probably  un- 
realistic. However,  the  initial  courses 
and  speeds  can  be  user  inputs. 

(b)  The  internal  TL  model  may  lack  the  detail 
and  accuracy  needed  for  fluctuation  studies. 
Range  or  azimuth-dependent  environments 
(sound  speeds,  bathymetry)  cannot  be  modeled. 

(c)  The  ship  source-level  model  is  unrealistic. 
Since  the  ship  speeds  and  lengths  are  at 
present  available  in  the  program,  the  inclu- 
sion of  a source  function  dependent  on  these 
parameters  could  be  accomplished  easily. 

(d)  There  is  no  mechanism  for  perfoi'ming  multiple 
replications  and  no  statistical  analysis 
routine.  But  again,  the  addition  of  such  a 
package  would  be  stral ght foi'ward . 
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(e)  The  TL,  receiver  response,  source  levels, 

ship  speeds  and  courses  are  time-independent. 

The  cost  ot  running  NABTAM  and  the  other  Brute-Force  models 
is  also  an  important  consideration,  and  there  are  only  rough 
guidelines  on  the  number  of  replications  required.  On  the 
other  hand,  NABTAM  can  run  several  receiver  locations,  beam 
patterns,  etc.  at  one  time.  Finally,  this  and  other  Brute- 
Force  models  may  be  the  most  efficient  means  for  determining 
short-term  statistics,  level  crossing  properties,  etc. 
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2.6  DISCRETE  SHIPPING  BEAM-NOISE  MODEL  (DSBN) 

2.6.1  Background 

Name:  (Discrete  Shipping  Beam-Noise  Model) (DSBN) 

Developer:  Science  Applications,  Inc., 

C.  W.  Spofford,  R.  G.  Stieglitz, 

H.  M.  Garon,  R.  C.  Cavanagh 

Sponsor:  Office  of  Naval  Research,  Code  431 

Previous  Applications:  Studies  for  APL/ Johns 

Hopkins  University,  LRAPP, 
Institute  for  Defense 
Analysis,  and  ONR  431. 

Published  Documentation:  Reference  F-1 

2.6.2  General  Approach 

Just  as  for  NABTAM,  the  DSBN  model  is  a Brute- 
Force  model  based  on  equation  (1-1).  The  beam-noise  time 
series  is  generated  from  component  submodels  for  the  sur- 
face ships,  TL,  and  receiver.  The  attributes  which  dis- 
tinguish DSBN  from  most  of  the  other  Brute-Force  models 
are : 


(a)  the  surface  ship  model,  which  uses  "Poisson 
lanes"  like  those  of  the  BTL  Analytic  model 
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(b)  a receiver  model  which  can  account  for  the 
vertical  as  well  as  horizontal  ship  signal 
arrivals 

(c)  a computer  structure  designed  for  multiple 
replications  and  multiple  beam  patterns 

(d)  an  extensive  statistical  analysis  package 

The  model  was  developed  for  investigating  the  fluctuations 
caused  by  sources  moving  through  the  TL  interference  pat- 
terns and  beam- response  azimuths,  so  that  attention  is 
paid  to  the  details  of  the  TL  and  the  beamformer. 

The  model  is  used  to  generate  Monte-Carlo  simu- 
ulations  of  noise  time  series  and  is  quite  general  in  its 
assumptions.  However,  TL  and  shipping  lanes  must  be  con- 
structed manually  and  are  considered  program  inputs.  The 
analysis  package  can  yield  the  usual  first  and  second 
order  (time)  statistics,  as  well  as  level-crossing  pro- 
perties, spectra  of  the  fluctuations,  ensemble  statistics 
over  replications,  beam-to-beam  cross-correlation  functions 
and  spectra. 

DSBN  is  programmed  in  FORTRAN  IV  and  is  presently 
run  on  CDC-6000-series  computers.  Cost  of  a 12-hour  repli- 
cation with  9 beam  patterns  and  200  ships  is  under  $5. 
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Comprehensive  statistical  analysis  of  ten  such  replications 
costs  about  $20. 


2.6.3  Model  for  Ships 

The  module  which  constructs  surface  ship  posi- 
tions and  parameters  can  be  viewed  as  having  two  parts. 

In  the  first  part,  the  source  levels,  courses,  and  speeds 
are  generated  as  realizations  of  random  variables.  The 
second  part  simply  tracks  each  ship's  position  and  com- 
putes bearing  and  range  to  the  array  as  functions  of 
time.  All  sailing  is  "rectangular,"  i.e.  on  a flat  earth 
and  not  on  great-circle  routes. 

The  determination  of  the  ships  and  their  para- 
meters proceeds  as  follows.  The  user  supplies  constraints: 

• Shipping  "lanes"  are  specified  by  the  dis- 
tributions of  speed  and  course  and  initial 
position,  as  well  as  an  inter-arrival  time 
interval  (expected  time  between  arrivals 
of  ships  across  a line  perpendicular  to 
the  mean  lane  course).  Lanes  are  of  finite 
length  and,  together  with  the  TL  function, 
define  the  boundaries  of  the  basin. 

• Source  levels  depend  on  the  random  speed, 
but  also  on  another  random  variable  (thought 
of  as  length)  whose  distribution  is  an  in- 
put. On  a given  lane  all  source-levels 
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are  drawn  independently  from  the  same  dis- 
tributions. The  source  intensity  does  not 
depend  on  vertical  or  azimuthal  transmission 
angle.  Note  that  the  source-level  model 
is  the  same  as  SIAM'S. 

The  present  version  of  the  DSBN  Model  uses  Poisson-distri- 
buted arrival  times,  so  that  the  lane  is  expected  to  have 
an  approximately  uniform  distribution  of  ships  on  the 
lane  (i.e.,  constant  density).  The  program  uses  a random 
number  generator  and  the  distribution  functions  to  pro- 
duce a single  replication  of  a shipping  field  with  speeds, 
courses,  locations,  and  source  levels  for  each  ship. 

Once  selected,  the  course,  speed  and  source  level 
of  a ship  remain  constant  for  the  duration  of  the  replica- 
tion time  (e.g. , 10  hours).  For  each  subsequent  replica- 
tion, the  selection  process  is  repeated. 

Note  that  the  shipping  model  described  here  is 
the  same  as  that  of  the  BTL  Analytic  model,  and  in  a sense 
consistent  with  that  of  the  WAGNER  model.  It  can  also  be 
shown  that  the  time- Independent  ship  locations  of  the  USl 
DX>del  are  consistent  with  this  Poisson-urrlval  approach. 
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2.6.4  Model  for  TL 

DSBN  requires  as  input  the  deterministic  TL  as 
a function  of  range  from  the  receiver  for  the  proper 
source/ receiver  geometries  and  frequency,  and  for  ranges 
to  the  limit  of  the  ocean  basin  (the  greatest  distance 
to  a target  or  surface  ship  source).  The  present  con- 
figuration of  DSBN  allows  for  the  use  of  the  vertical 
arrival  structure  (i.e.,  the  loss  as  a function  of  range 
and  vertical  angle  at  the  receiver)  in  order  to  model  the 
three-dimensional  response  of  an  array.  The  model  can 
also  treat  TL  as  a function  of  bearing;  but  for  more 
than  a few  different  angles,  running  time  and  computer 
storage  can  be  large.  DSBN  cannot  utilize  a TL  function 
which  depends  on  time.  The  incorporation  of  such  a fea- 
ture amounts  to  a bookkeeping  problem  and  could  be  han- 
dled if  the  detail  and  computer  expense  were  warranted. 

There  is  no  real  limit  to  the  amount  of  detail 
permitted  in  the  TL  input.  For  studies  of  source-motion- 
induced  fluctuations  caused  by  multipath  interference, 
the  range  resolution  in  the  TL  model  depends  on  acoustic 
frequency,  processor  integration  time,  and  the  speed  of 
the  sources  (i.e.,  the  velocity  component  radial  to  the 
receiver).  In  a special  study  reported  in  Ueference  F-1, 
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the  TL  sampling  rate  appropriate  to  25  Hz,  2-minute  inte- 
gration times,  and  15-20  knot  speeds  was  found  to  be 
about  0.2  miles.  The  maximum  range  for  ship  contributions 
in  the  basin  was  assumed  to  be  500  miles.  Thus,  a typical 
TL  table  for  low  acoustic  frequencies  and  open  ocean  has 
500  X 5 = 2500  entries  for  each  source/ receiver  depth 
combination.  It  should  be  clear  that  to  incorporate 
time  or  bearing-dependent  TL  requires  tables  with  perhaps 
2500  X 10  to  2500  x 1000  entries.  If  arrival  structure 
were  included,  multiply  these  figures  by  10-100. 

2.6.5  Receiver  Model 

The  receiver  location  is  fixed,  and  does  not 
change  with  time.  For  the  present  configuration  of  the 
DSBN  Model,  the  receiver  module  is  simply  a functional 
giving  the  Intensity  response  of  the  array  to  plane  wave 
arrivals  (l.e.,  beam  pattern).  Since  the  usual  pi'oblero 
deals  with  a horizontal  line  array,  the  response  is  often 
given  in  terms  of  azimuthal  arrival  angles.  However, 
to  investigate  the  effects  of  the  horizontal  array's 
vertical  directivity  away  from  broadside,  the  simulatU>n 
model  can  accommodate  a response  function  which  depends 
on  both  azimuthal  and  vertical  arrival  angles. 
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The  contributions  at  the  beamformer  output  of 
the  various  ships  contributing  to  the  noise  field  are 
added  on  a random  phase  (incoherent)  basis.  All  temporal 
processing,  filtering,  etc.  are  handled  implicitly  in 
the  source  level  and  TL  functions  for  the  basic  time 
samp  I e . 


For  computation  purposes,  the  user  specifies 
the  fixed  location  and  depth  of  the  array  (for  the  geo- 
metry of  the  sources  and  the  TL)  and  provides  up  to  ten 
array  response  functions  corresponding  to  different 
steering  angle  orientations,  shading,  physical  deforma- 
tions, or  whatever.  The  DSBN  Model  simulates  beamformer 
output  for  each  response  function  or  beam  pattern  by 
modifying  the  intensity  arrivals  from  ship  sources 
accordingly.  The  code  also  records  the  number  of  sources 
on  the  "niain  beam",  defined  at  input. 

For  ease  of  operation,  two  idealized  array  re- 
sponse functions  are  programmed  as  optional  subixiutines 
for  DSBN: 


(a)  Shaded  Horizontal  Line  Array  - Azimuthal 
Uesponse  Only.  The  function  approximates 
the  response  of  a shaded  horizontal  line 
array  (uniform  spacing)  with  prescribed 
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main  beamwidth  and  sidelobe  suppression. 

The  sldelobes  have  structure  approximating 
Dolph-Chebyshev  shading. 

(b)  Shaded  Horizontal  Line  Array  - Vertical 

and  Azimuthal  Response.  As  in  (a),  except 
that  sidelobes  are  completely  suppressed. 

Given  that  undistorted  plane-wave  response  suf- 
fices, the  primary  limitation  of  the  program  is  that  the 
response  function  and  array  location  cannot  change  with 
time.  Such  features  can  be  added  with  minimal  modifica- 
tions to  simulate,  for  example,  the  response  of  a transiting 
towed  array  which  is  chiinging  orientation  and  suffering 
from  physical  deformation  (wiggles). 

A more  basic  limitation  of  this  and  other  mtidels 
reviewed  here  Is  that  the  array  model  does  not  apply  directly 
to  a predicted  acoustic  field  which  has  not  been  decixiposed 
into  plane  waves  (e.g.,  output  of  the  PK  model).  The  most 
efficient  way  to  deal  with  this  is  to  include  the  array 
response  in  the  transmission  loss  by.  in  effect,  immersing 
the  array  elements  into  the  field  and  computing  the  beam- 
former  algorithm  directly.  This  has  been  done  for  a single 
source  and  could  certainly  be  extended  to  the  noise  case. 
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2.6.6  Details  of  the  Calculation 

Once  the  shipping  field,  TL,  and  array  response 
are  determined,  the  noise  is  computed  at  discrete  time 
steps  for  each  of  ten  (or  fewer)  array  response  functions 
according  to  equation  (1-1).  Readings  are  taken  at  pre- 
specified time  points  (e.g.,  every  minute)  for  a selectable 
time  span  (e.g.,  30  hours).  The  ship  model  provides  the 
ship  positions  and  source  levels;  the  environment  and  TL 
models  provide  the  transmission  loss  from  each  source  to 
i,he  receiver,  and  the  array  model  provides  the  beam  response. 
The  modular  form  and  flow  of  the  computer  routine  are 
sketched  in  Figures  F-1  and  F-2.  Note  that  only  one  re- 
ceiver depth  and  frequency  is  treated  per  replication. 

The  statistical  analysis  packages  operate  on 
tlw  simulated  time  series  and  are  described  next. 


2.6.7  Output  and  Analysis 

For  each  replication,  the  output  of  the  DSUN 
model  consists  of  beam-noise  time  series  for  up  to  ten 
beam  patterns  ("beams'*)  and  for  a single  frequency. 

The  results  are  then  analysed  by  several  computer  packages  - 
which  operate  on  single  or  multiple  replications. 
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2 . 6 . 7 . 1 Smt  1st, leal  Analysis  Routine 

A ttenerul  statistical  packURo  has  been  constructed 
to  study  array  time-series  data.  It  operates  on  the  matrix 

fj.  t^) 

where  f ^ , and  tj^  are  interpreted  as  the  discrete  beam 
pattern  indices,  frequencies,  and  time  steps,  respect i\rely. 

N is  the  noivio  level  (in  dB).  The  fo.UowinR  calculations 
art-  performed. 

(u)  Histograms  are  constructed  and  plottt'd 
for  any  raiiRe  of  i,  .i  and  k xo  specified 
resolution.  Likewise,  the  mean,  variance, 
skewness,  kurtt>sis,  and  decile.s  are  fi'und. 

(b)  For  two  of  the  three  independent  variables 
fixed . 

• The  series  is  plotted 

• A "stationnrity"  test  is  periorn^ed  by 
divldiuR  the  series  into  any  number 
of  equal  parts  and  then  applying  (u) 
to  each  part 

• The  sainjjle  autocovnriance  function  is 
eomiJUted  and  plotted 

• The  autocorrelation  function  is  esti- 
iwited 
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• An  FFT  is  applied  to  the  autocorrelation 
function  to  estimate  the  power  spectral 
density 

(c)  For  one  independent  variable  fixed,  the  two- 
dimensional  autocorrelation  function  is  esti- 
mated and  output  in  matrix  form. 

(d)  For  separation  (lag)  in  one  variable,  en- 
sembling  over  the  second  and  for  the  third 
fixed,  the  cross-correlation  function  is 
found. 

(e)  For  one  variable  fixed,  one  separated,  and 
one  lugged,  the  joint  density  function  for 
the  separated  variables  is  estimated.  The 
histogram  Is  found  and  multivariate  moments 
calculated. 

(f)  A Lllllefors  Test  (see,  e.g. , Hef.  F-2)  for 
got>dness-of-f i t can  be  applied  to  the  sample 
histogram  to  find  best  Gaussian  fit  and  test 
at  confidence  levels  of  0.95  and  0.99. 

(g)  The  logarithmic  transfoj’raation  of  the  Log- 
Normal,  Non-Central  Chi-Square,  and  Chi- 
Square  distributions  are  tested  against  the 
sannle  distribution  (for  best  estimate  of 
parameters  based  on  median  and  one  or  more 
other  percentile  vtoints)  at  levels  0.95 
and  0.99  with  the  Kolmogorov  Test  for  fit 
(see  Uef.  F-2).  Graph.s  of  the  sample  and 
fitted  functions  are  plotted. 
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(h)  A simplified  test  for  ergodicity  calculates 
ensemble  statistics  in  two  directions  (e.g. , 
in  t and  then  in  replicas  for  f and  i}*  fixed) 
and  compares  sample  distribution  functions 
at  the  0.95  and  0.99  levels  with  the  Smirnov 
Test  (Ref.  F-2). 


2 . 6 . 7 . 2 Detector  (Level-Crossings)  Analysis  Routine 

A computer  package  has  been  designed  to  model 
several  types  of  detectors.  Input  consists  of  a time  series 
plus  relevant  parameters.  The  input  time  series  is  converted 
to  a time  history  of  detect /no  detect  states  according  to 
the  following  algorithms. 

Continuous  Threshold  detector 

Given  a time  series  X(t)  . a threshold  TH, 
and  a time  interval  T (holding  time),  score 
a detection  at  time  t^^  if  X(l)  TU  con- 
tiuut'usly  for  t^  - T ^ t ^ t^^. 

(b>  Union  of  Continuous  Threshold  detector 

This  is  a generalization  of  (a)  except  that 
a sequence  of  thresholds  and  associated  holding 
periods  tTH^  Tj)  comprise.s  the  input.  Then 
detection  occurs  if  the  signal  X(t)  has  con- 
tinuously exceeded  a given  threshold  for  the 
associated  holding  period  for  any  menjber  of 
the  sequence  (THj . 
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(c)  Intensity  Average  detector 

Given  a time  series  X(t)  , a threshold  TH, 
and  an  averaging  time  T,  construct 


where  the  sum  extends  over  all  times  tj  such 
that  tQ  - T < tj  £ tQ.  Score  a detect  at 
time  Iq  if  10  logjQ  S(t)  ^ TH,  i.e. , if  the 
intensity-a 'eraged  X(t)  exceeds  TH. 


Union  of  Intensity  Average  detector 


This  generalizes  (c).  The  input  consists 
of  a sequence  of  thresholds  and  associated 
averaging  Intervals  (TH^,  T^) . A detection 
occurs  at  time  t^  if  the  intensity  average 
over  any  one  of  the  averaging  periods,  T^ , 
exceeds  the  associated  threshold.  TH^ . 


(e)  N Out  of  M detector 

Given  a time  history  record  X(t)  . a 
threshold  TH,  and  integers  N and  M.  a 
detection  occurs  if  X{t)  has  exceeded 
TH  for  at  least  N out  of  M time  points 
immediately  preceding  and  including  t^j. 


» • 

f Any  one  of  these  detectors  yields  a lime  history  of  detect/ 

no  detect.  Various  statistics  are  then  calculated  and 
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displayed,  Including  the  distribution  of  detect  (holding) 
times,  no-detect  times,  associated  moments,  order  statis- 
tics, waiting  times  to  cross  a threshold,  a complete  history 
of  beam-free  and  beam-occupied  times  and  associated  moments. 


2.6.8  Computer  Implementation 

The  DSBN  model  and  associated  analysis  packages 
are  programnied  in  FORTRAN  IV  and  have  been  run  on  CDC  6000- 
series  computers.  The  shipping  lanes,  TL,  and  beam  patterns 
(if  different  from  the  two  listed  in  Subsection  2,6.5)  are 
program  Inputs  and  are  not  automatically  generated. 

Costs  depend  on  the  number  of  replications,  ships, 
beam  patterns,  time  period  and  sampling  rate,  TL  detail,  etc., 
and  of  course  on  the  analysis.  Examples; 


• 200  ships.  12  hours.  2 minute  samples. 

0 beam  patterns,  1 frequency.  2 nm  average 
TL.  horizontal  reupiinse  only:  $5/replicut ion 

• The  same  case  for  lO-mlnute  sampling; 

$l/repl ication 

• Analysis  of  8 12-hour  replications  such  as 
given  above  (see  the  Table  of  Section  3) 

$20. 
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Core  for  such  cases  is  typically  under  lOOK  (octal)  or  32K 
(decimal)  words. 

2.6.9  Evaluation 

Outputs  of  the  DSBN  model  are  at  present  being 
compared  with  Square-Deal  and  other  measurement  data  for 
which  there  are  shipping  and  environmental  data.  No  con- 
clusion can  be  made  at  this  time. 

2.6.10  Significant  Advantages  and  Disadvantages 

Since  DSBN  was  designed  to  investigate  detailed 

beam-noise  fluctuations  caused  by  source  movements,  it  has 
such  special  features  us;  a computer  structure  which  accomo- 
dates multiple  beam  patterns  and  replications,  a receiver 
model  that  cun  simulate  multipath  beam  splitting,  an  exten- 
sive time  series  analysis  package  to  generate  first  and 
second  order  distribution  functions  and  »n\->«u>nts  as  well  as 
level  crossing  statistics  and  fluctuation  spectra,  with  time 
and  ensemble  averaging. 

The  model  has  a number  of  intrinsic  limitations. 
The  Poisson  ship-arrival  model  has  not  been  verified;  'he 
Tl.,  source  level,  etc.  are  t i me- Independent ; it  is  core- 
consuming  to  have  more  than  a few  different  TL  curves  as 
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functions  of  bearing;  only  one  frequency  and  one  receiver 
depth  are  treated  in  a single  run;  the  beam  response  functions 
are  limited  to  special  form.  The  modular  form  of  DSBN  allows 
for  relatively  easy  elimination  of  these  deficiencies. 

The  cost  of  running  DSBN,  and  the  other  Brute- 
Force  models  as  well,  is  an  important  consideration.  More- 
over, there  are  at  present  only  rough  guidelines  on  the 
number  of  replications  required  to  obtain  a meaningful  sta- 
tistical sample.  On  the  other  hand  we  note  that  once  the 
shipping  lanes  for  a basin  are  established,  the  model  can 
be  run  at  multiple  locations  without  additional  preparation 
of  ship  data.  Furthonnore,  for  determining  the  short-term 
statistics,  level-crossing  properties,  etc.,  the  Brute-Force 
approach  may  be  the  most  efficient  one. 
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BEAMPL 

Background 

Name:  (BEAMPL)  ("PL"  stands  for  "Program  Library") 

Developer:  Office  of  Naval  Research  (AESD), 

C.  W.  Spofford,  R.  G.  Stieglitz, 

H.  M.  Garon,  R.  C.  Cavanagh 

Sponsor:  Office  of  Naval  Research  (LRAPP) 

Previous  Applications;  Studies  for  LRAPP,  PM-4 

Published  Documentation:  None 

Current  Residence:  NORDA,  Code  31)0 


t 

I 

i 


2.7.2  General  Approach  and  Sunmmry 

The  BEAMPL  noise  nK>dt>l  is  nearly  identical  to 
DSBN.  In  fact,  the  basic  ideas  of  BEAMPL  were  modified 
and  extended  in  the  construction  of  DSBN,  so  that  wt' 
limit  this  description  to  an  identification  of  the  parts 
of  BEAMI>L  which  differ  from  DSBN. 


I 


• Ship  source  intensities  are  deleministic 
constants  for  each  class  of  ship. 


1 

I 


• Ship  lanes  and  movements  are  as  in  DSBN 
except  that  all  ships  in  a lane  have  the 
same  deterministic  constant  speed. 
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Tl.  is  Independent  of  azimuth  and  the  effects 
of  vertical  arrival  structure  are  not  in- 
cluded. 

The  receiver  module  differs  from  that  of  DSBN 
in  that  it  allows  only  one  beam  pattern,  which 
depends  only  on  the  azimuthal  arrival  angles 
and  has  "spotlight”  response  (the  response 
is  unity  on  the  prespecified  main  beam,  and 
zero  elsewhere). 

The  output  is  the  same  as  for  DSBN,  but  the 
analysis  packages  are  limited  (see  the  Tables 
in  Section  3). 

The  computer  routine  for  DEAMPL  is  sub- 
stantially different  from  that  of  DSBN  and 
is  somewhat  less  efficient. 
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2.8  SIAM  I 

2.8.1  Background 

Name:  (Simulated  Ambient  Noise  I)  (SIAM  I) 

Developer:  Naval  Research  Laboratory, 

Samuel  V.  Marshall,  John  J.  Corny n 

Sponsor:  Office  of  Naval  Research, 

LRAPP 

Previous  Applications:  Ex'aluation  with  lOMEDEX 

and  other  data  for  LRAPP 

Published  Documentation;  References  H-1  and  K-2 

2.8.2  General  Approach 

There  are  two  models  which  share  the  name  '‘SIAM," 
In  this  section  we  consider  the  earlier  version,  and  call  it 
"SIAM  I"  or  just  "SIAM."  The  model  is  not  actively  used  at 
the  present  time,  but  It  has  features  worth  reviewing.  In 
the  next  section,  the  successor,  called  "SIAM  11."  will  be 
described,  it  is  substantially  different  from  SIAM  1. 

*The  name  "SIAM"  originally  referred  to  a driver  for  the 
NRL  model  which  randomly  initializes  the  shipping  param- 
eters. controls  the  Monte-Carlo  simulations,  av  , calcu- 
lates the  ensemble  statistics.  Ic  have  used  "SIAM"  here 
as  the  name  for  the  entire  noise  routine. 
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SIAM  is  a Brute-Force  model  designed  to  predict 
ship-generated  noise  over  the  band  20-120  Uz.  It  is  based 
on  Equation  (1-1).  The  team-noise*  time  series  are  calcu- 
lated with  ctwiponent  submodels  for  the  ships.  TL,  and 
receiver.  Notable  is  the  ship  scenario  in  which  ships  travel 
on  specified  courses,  but  reflect  from  the  basin  boundary 
and  return  (a.s  in  the  WAGNER  Bounded  noise  process).  Initial 
ship  courses,  speeds,  source  levels  are  input  or  selected 
from  distributions,  but  remain  constant  thereafter.  It  isi 
of  special  interest  that  the  TL  has  a random  component  (as 
in  the  USl  model)  which  is  chosen  independently  for  each 
source  and  each  time  step  from  a common  distribution.  The 
model  iKDploys  great-circle  geometry  and  can  predict  the 
simuUantH>us  output  of  several  different  deterministic  and 
t Ime-i  Independent  array-response  functions.  Standard  oper- 
ating procedure  for  SIAM  is  to  generate  many  replications, 
so  that  enstHttble  statistics  are  emphasized. 

The  model  has  been  run  on  the  CDC  3800  computer 
at  NRL  and  is  programmed  in  FORTRAN  IV. 


As  noted  below,  only  idoalizcMl  "SrOTLlCiHT'’  beam  patterns 
are  used. 
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2.8.3 


Model  for  Ships 


A finite  basin  is  defined  in  advance,  using  a 
maximum  of  100  points  to  define  the  boundary  vertices.  Then 
ship  locations,  courses  and  speeds  are  initialized  by  class 
either  deterministically  or  from  random  distributions.  Once 
these  parameters  are  established,  the  ships  travel  on  great- 
circle  paths  at  constant  speed  and  with  constant  source 
level.  When  a ship  encounters  the  basin  boundary,  it  is 
reflected  specularly. 


The  source  intensity  is  a random  variable  for 
each  ship  and  is  determined  from  a distribution  appropriate 
to  its  class.  It  depends  on  class,  frequency,  speed,  and 
length,  the  latter  two  of  which  are  uniform  or  normal  vari- 
ables. Default  values  for  source  level  at  50  Hz  are 


10  log  (SL)  = SLq  + 60  log  (V/V)  + 20  log  (L/L)  dB 
where  SLq  Is  normal  with  mean  160  and  o = 3 (dB) 

V =*  16,  V normal  with  mean  16  and  a = 5 (knots) 

L = 525,  L normal  with  mean  525  and  a = 170  (feet) 


Source  levels  are  Independent  of  aspect,  vertical  angle,  and 
time.  A maximum  of  250  ships  is  allowed. 
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2.8.4 


Model  for  TI, 


SIAM  requires  as  Input  a <ietermlnisiic , time- 


independent  TL  as  a function  of  range,  a'/imuth,  frequency, 


and  receiver  depth.  In  lieu  of  a user-specified  table,  the 


program  can  generate  TL  of  form  A+BlogU.  The  model  uses 


a special:  "logarithmic"  routine  for  interpolating  TL  in 


range  and  azimuth.  Moreover,  all  geometry  is  spherical. 


The  only  limit  in  the  amount  of  detail  permitted  for  the  TL 


input  is  one  of  computer  core  and  running  time. 


Current  restrictions  are: 


• a maximvan  of  10  TL-versus-range  ourvt's 
for  each  frequency  and  lu'ceiver  depth 


• a maximum  of  570  points  for  each  TL 


curve 


• a maximum  of  30  frequencies,  29  depths 


each 


SIAM  also  adds  a random  fUictuation  component 


to  the  TL.  The  values  are  select  <‘d  indt'pendt'ni  ly  at  each 


time  step  fiann  an  input  distribution  idt*fault  is  a normal 


distribution  with  mean  zero  and  o = 1.5  dB). 
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2.8.5  Receiver  Model 
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The  receiver  location  is  an  input  which  does  not 
change  with  time.  The  model  is  designed  to  give  the  response 
of  an  idealized  horizontal  array,  but  does  not  account  for 
vertical  arrival  structure  for  off  broadside  steering  angles. 
The  array  response  is  treated  as  a "spotlight"  beam  pattern: 
with  unit  intensity  responses  on  the  "main  beam"  and  zero 
elsewhere.  Tae  model  assumes  that  all  energy  arrives  as 
horizontal  plane  waves,  perfectly  coherent  and  undistorted, 
from  the  source's  azimuth.  Contributions  from  the  various 
ship  sources  are  summed  on  a random  phase  basis  (incoherently) 
and  the  effect  of  temporal  processing,  filtering,  etc.  on  a 
basic  time  sample  are  included  in  the  source  level  and  TL 
functions. 


SIAM  is  structured  to  yield  simulations  of  the 
simultaneous  noise  output  for  multiple  beams.  Hence,  the 
input  consists  of  several  beam  patterns,  i.e.,  several  main 
beams,  but  no  side lobes. 

2.8.6  Details  of  the  Calculation 

Once  the  ships,  TL,  and  array  response  are  estab- 
lished, as  discussed  above,  the  noise  is  conjputed  at  discrete 
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time  steps  according  to  equation  (1-1).  Readings  are  taken 
at  specified  time  points  for  a selectable  time  span  but  with 
a maximum  of  4000  samples  per  case.  The  computer  flow  is 
summarized  in  Figure  H-1. 

2.8.7  Output  and  Analysis 

The  basic  SIAM  output  is  the  time  series  of  noise 
as  a function  of  beam,  frequency,  and  depth.  The  usual  oper- 
ating procedure  is  to  perform  many  replications  over  the 
various  ^luctaators  (TL,  SL,  ship  parameters).  A special 
analysis  routine  prints  or  plots  ship  locations  and  noise 
time  ae''ies,  and  calculates  the  "grand  ensemble"  moments  of 
the  noise,  l.e,,  averages  are  over  time  and  replication,  as 
well  as  the  sample  one-dimoasional  density  function  (all  for 
noise  levels  in  dBs). 

The  program  y..elds  complete  histories  of  ships 
(including  plots)  as  well  as  percentages  of  beam-free  and 
beam-occupied  tines.  A valuable  feature  of  the  SIAM  output 
is  the  identification  of  the  ship  which  r.akos  the  greatest 
contribution  at  each  time  step  and  a list  of  ships  in  order 
of  thair  importance  ov-^r  a specified  time  interval. 
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Guidelines,  based  on  experimental  runs  of  SIAM 
for  canonical  cases,  have  been  established  for  the  number 
of  Monte-Carlo  replications  required  to  obtain  "convergence." 

2.8.8  Computer  Implementation 

SIAM  is  programmed  in  FORTRAN  IV  and  has  been 
operated  on  the  NRL  CDC  3800  computer.  The  shipping  parsun- 
eters,  TL,  beam  patterns  are  program  inputs  and  are  not  auto- 
matically generated. 

Costs  depend  on  the  number  of  replications,  fre- 
quencies, depths,  beams,  ships,  sample  times,  etc.  The  author 
reports  that  a typical  run,  as  reported  in  Ref.  H-2,  would 
require  no  more  than  5 or  10  minutes  on  the  CDC  3800  computer. 

2.8.9  Evaluation 

SIAM  model  output  has  been  compared  with  beam- 
noise  data  from  controlled  experiments  (see  Ref.  H-2).  The 
uncertainties  in  the  measurement  data,  input  data,  and  analy- 
sis preclude  a thorough  evaluation  of  the  model  but  agree- 
ment between  measured  and  predicted  noise  distribution  func- 
tions was  generally  good. 
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2.8.10  Significant  Advantages  and  Disadvantages 

SIAM  I was  designed  to  generate  Monte-Carlo  esti- 
mates of  one-dimensional  distributions  of  noise.  Its  treat- 
ment of  ship  movement  and  its  computer  structure  are  tailored 
to  this  application.  Noteworthy  features  are  the  great-circle 
ship  sailing  and  the  additive  TL  fluctuation  term. 

The  model  has  certain  Intrinsic  limitations.  The 
deterministic  part  of  the  TL,  the  receiver  location  and 
response,  the  source  level,  ship  speeds  and  courses  are  all 
independent  of  time;  the  vertical  arrival  structure  at  the 
array  is  not  accounted  for;  computer  storage  limits  the 
number  of  ships;  the  ship  reflection  model  may  be  unrealistic; 
the  analysis  package  is  limited  Furthermore,  the  addition 
of  a TL  fluctuation  terra  independently  at  each  time  step 
precludes  the  possibility  of  calculating  realistic  temporal 
autocorrelation  or  higher  order  moments.  Each  of  these  could 
be  eliminated  in  a straightforward  way. 

The  cost  of  running  SIAM,  and  the  other  Brute- 
Force  models  a.s  well,  is  an  important  consideration.  How- 
ever, the  model  can  be  run  for  several  frequencies  and 
receiver  depths  and  beams  at  one  time,  and  once  the  ship 
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field  is  input,  it  can  be  used  for  other  receiver  sites  in 
the  basin.  Finally,  for  determining  short-term  statistics, 
level  crossing  properties,  etc.,  the  Brute-Force  approach 
may  be  the  most  efficient  one. 
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2.9  SIAM  II 

2.9.1  Background 

Name;  (Simulated  Ambient  Noise  II)  (SIAM  II) 

Developer:  Naval  Research  Laboratory, 

S.  C.  Wales.  S.  W.  Marshall 

Sponsor:  Office  of  Naval  Research,  LRAPP 

Previous  Applications:  Evaluation  with  EPAC 

and  other  data  for  LRAPP 
(References  A-1,  I-l  and 
1-2) 

Published  Documentalon:  In  preparation  (User's 

Manual/JUA  Article). 

See  Reference  1-2. 


2.9.2  General  Approach 

SIAM  II  is  the  successor  to  SIAM  I (Section  2.8). 
It  is  a Brute-Korce  model,  designed  to  provide  many  replica- 
tions of  surface-ship  noise  for  horizontal  array  syattmis. 
especially  narrow-beam  systems.  It  differs  frtan  the  other 
noise  models  in  its  approach  to  simulating  ship  traffic, 
but  does  calculate  beam  noise  from  equation  (1-1). 
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The  unique  aspect  of  the  SIAM  II  approach  is  that 
it  calculates  many  (32)  sample  time-series  in  one  run,  using 
a special  ship-tracking  routine  which,  in  effect,  cycles  ships 
through  the  beams  within  annular  regions.  This  is  done  to 
solve  the  problem  of  creating  new  ships  as  others  leave  the 
basin  and  to  provide  multiple  samples  from  the  same  batlc 
ship  scenario.  Source  levels  and  speeds  are  con.^tant  in  time, 
but  initialized  from  random  populations. 


The  TL  is  an  input  table  for  each  azimuthal  sector, 
but  does  not  depend  on  vertical  angle.  It  has  a fluctuation 
component,  uncorrelated  in  time.  The  array  response  is  modeled 
with  a single  deterministic  beam  pattern. 


The  model  is  implemented  on  a PDP  11-45  computer 
and  is  prograimaed  in  FORTRAN.  The  present  configuration  makes 
for  lengthy  running  times  and  tedious  data  entry. 


2.9.3 


Model  for  Ships 


The  SIAM  II  ship  traffic  model  is  different  from 
all  of  the  others  reviewed  here,  and  we  spend  some  time 
describing  it.  An  ocean  "basin,"  centered  at  the  fixed 
receiver  location,  is  divided  into  sectors  and  subsectors 
as  shown  below. 
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A woctor 


A Subsector 


The  sector  aperture  is  determined  from  the  array  beam  pattern 
<e.g.,  corresponding  to  the  main  beam  width  and  perhaps 
principal  sidelobes),  while  the  subsector  extent  is  chosen 
to  be  consistent  with  shipping  lane  widths  (usually  50-200  nm) 
and  the  basin  environment.  A maximum  of  16  sectors  is  allowed, 
but  the  number  of  subsectors  is  unlimited.  All  geometry  assumes 
a flat  earth.  As  discussed  below,  the  contribution  to  the 
noise  field  from  each  of  the  subsectors  is  calculated  sequen- 
tially. 
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For  a given  subsector,  the  computer  routine  then, 
in  effect,  constructs  vp  to  3S  copies  of  tb^  subsector  and 
internally  generates  an  annulus  from  the  copies,  centered  at 
the  receiver  location; 


Original 

Subsector 


Ships  are  Initially  distributed  uniformly  in  the  annulus 
according  to  the  input  subsector  density.  Ship  speeds  (and 
lengths)  are  chosen  independently  from  a nearly-normal  distr. 
butlon  for  prespecified  mean  and  variance.  Ship  courses 
are  drawn  from  a uniform  distribution  over  {0^ , 360^').  Ship 
source  levels  may  be  inputs  or  may  follow  a speed/ length 
rule  as  in  DSBN  or  SIAH  1.  Once  selected,  these  ship  param- 
liters  remain  constant  for  the  time  span  of  one  replivalioo 
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(up  to  128  tiine  steps).  The  ships  are  then  tracked  in  time 
around  the  annulus,  following  straight-line  paths  and  reflec- 
ting specularly  from  the  boundaries  as  encountered: 


The  ship  history  within  each  subsector  copy  is  treated  as  u 
sample  o^  ship  traffic  for  the  original  subsector,  so  that 
32  samples  are  obtained  at  once.  This  is  the  objective  of 
the  approach.  If  the  sector  aperture  exceeds  11. 25^,  i.e,, 
360^/32.  then  fewer  than  32  samples  are  obtained  in  a single 
run  and  more  runs  are  made  until  32  samples  are  obtained. 
Separating  the  sector  copies  is  allowed,  but  overlapping  them 
may  cause  undesirable  correlations  for  adjacent  copies. 


2-101 


I > SIAN  II 
Model 


The  above  procedure  is  repeated  for  each  subsector 
of  interest.  The  model  usually  concentrates  on  the  main  lobes 
of  an  array,  so  that  only  a few  sectors  need  be  treated  in  a 
single  model  run.  Full  360°  response  can,  however,  be  obtained. 


As  is  the  case  for  most  of  the  other  noise  models, 
source  intensity  is  independent  of  time  and  the  transmission 
angle. 


2.9.4 


Model  for  TL 


SIAM  II  uses  an  input  table  of  TL  versus  range  for 
each  of  as  many  as  16  sectors,  a single  frequency,  and  one 
receiver  depth.  As  an  alternative,  the  model  will  generate 
an  A BlogR  table  internally  for  each  sector.  The  amount 
of  detail  in  the  TL  function  is  selectable,  but  should  be  con- 
sistent with  the  time-sampling  rate,  i.e.,  the  rale  at  which 
the  time  series  of  noise  is  generated  and  the  sonar-system 
averaging  time. 


As  an  added  feature,  a fluctuation  component  can 
be  added  to  the  input  TL.  It  is  chosen  independently  for 
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The  program  is  not  structured  to  account  for  the 
vertical  arrival  structure  of  transmission. 


2.9.5 


Receiver  Model 


The  receiver  location  and  depth  are  implicitly 
defined  by  the  TL  and  ship  field,  and  do  not  change  in  time. 

The  receiver  is  assumed  to  be  a horizontal  array  with  response 
to  plane  wave  arrivals  in  the  horizontal  plane  determined  from 
an  input  beam  pattern.  No  discrimination  in  vertical  angles 
is  allowed.  Only  one  beam  pattern  (frequency,  array  depth) 
is  permitted  in  a single  replic  .tion.  As  an  option,  a special 
"spotligfit"  mode  employs  a single  "perfect"  beam  and  shortens 
computation  time. 

As  a ship  progresses  around  the  annulus,  its  contri- 
bution at  each  time  point  is  added  incoherently  to  the  appro- 
priate sector  copies.  Only  one  beam  pattern  is  allowed  for 
each  sector,  but  different  beam  patterns  may  be  input  for 
different  sectors  and  sidelobes  can  be  included  in  this  fashion. 
All  temporal  processing,  filtering,  etc.  are  treated  implicitly 
in  the  source  and  TL  functions  for  the  basic  tinie  samples. 
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tion. 

• The  user  inputs  a TL  table  for  each  sector, 
appropriate  for  the  single  frequency,  receiver 
depth,  and  surface-ship  source. 

• The  user  inputs  a beam  pattern,  time-step 
increment,  and  number  of  time  steps  (up  to 
128). 

• For  each  subsector,  the  user  inputs  the  ship 
parameters:  mean  number  of  ships,  mean  ship 
speed  and  length,  moan  (standard)  source  level. 


• The  program  begins  with  the  first  sector  and 
its  first  subsector 

• The  subsector  copies  (32)  and  annulus  are 
generated  Internally  (as  illustrated  above). 
The  program  initializes  ship  positions  in  the 
subsector,  ship  speed,  courses,  and  source 
levels  - using  the  appropriate  random-number 
generator. 
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Each  ship  then  moves  in  time  in  the  annulus 
for  128  time  points.  The  source  levels,  TL 
for  the  sector  and  range,  and  beam-response 
for  the  ship  bearing  are  combined  to  yield 
32  realizations  of  the  ship's  contribution  to 
the  noise  for  the  current  sector  as  a function 
of  time.  The  intensity  is  added  incoherently 
to  a cumulative  sum  of  noise  for  each  time 
point  and  the  first  sector.  These  values 
are  stored  in  a (copy  number  time)  array  of 
up  to  (32  X 128  points). 

Once  all  ships  for  the  current  subsector  have 
been  treated,  the  program  goes  on  to  the  next 
subsector,  and  continues  until  all  subsectors 
of  the  first  sector  are  completed.  The  result 
is  32  realizations  of  the  beam-noise  time 
series  for  ships  in  the  first  sector. 

Ship-noise  contributions  for  the  other  sectors 
are  then  calculated  in  sequence  in  the  same 
way  as  for  the  first  sector. 

As  each  sector  is  completed,  its  contributions 
to  the  beam  noise  (32  time  series)  are  added 
incoherently  to  the  cumulative  sum  stored  in 
the  32  X 128  point  array. 

Basic  output  consists  of  32  time  series  of  beam 
noise,  each  of  length  up  to  128  time  steps, 
as  well  as  ship  histories. 
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The  program  flow  is  summarized  in  Figure  1-1. 

Notice  that  in  the  innermost  loop  the  contribution  of  an 
individual  ship  is  determined  for  the  entire  time  span  and 
for  each  of  the  32  subsector  copies.  For  this  reason,  present 
program  size  limits  the  number  of  time  points  per  time-series 
sample  to  128.  Note  also  that  only  one  complete  beam  pattern 
is  allowed  per  computer  run. 


We  emphasize  that  new  ships  and  ship  tracks  are 
constructed  for  each  subsector's  contribution.  Hence,  there 
is  no  correlation  between  ships  of  one  subsector  and  ships  of 
another. 


Finally,  a constant  wind-noise  level  is  added  in- 
coherently to  the  resulting  ship-noise  time  series. 


2.9.7 


Output  and  Analysis 


The  basic  model  output  is  a set  of  32  beam-noise 


time  series  (or  several  sets  of  32  if  additonal  replications 
are  made).  An  analysis  package  operates  on  these  data  to 


yield: 


• Histograms  of  Noise  (1  dB  resolution) 


• Sample  Distribution  Functions  (in  dBs) 
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Set  up  Run  Files 

\ 

'-and  Transmission  Loss  Files 


Loop  in  Computations 


Loop  in  Sectors 


Define  Sector 
Input  TL(R) , 
Array  Response| 


Loop  in  Subsectors 


f ^Define  Subsector  , 
i Input  Ship  Density  , 
«.  and  Parameter 
Distributions 


Loop  in  Ships 


Select  Ship  Parameters, 
Calculate  Source  Level 


i Move  Ship  Within  Annulus 
j Calculate  Received  Level 
Loop  in  time  ! from  TL  and  Array  Response 
j for  appropriate  copies 

! and  add  to  time  series 


Add  Wind  Noise  to  all 
time  series 


Analyze  time  series 


Output  time  series  for 
further  analysis 


Figure  I-l.  SIAM  II  Program  Flow 
(Courtesy  of  S.  C.  Wales,  NRL) 
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• Mean  and  Variance  for  the  Noise  Levels  (dB) 

• Median  of  the  Noise  Distribution 

Temporal  autocorrelation  is  not  available  as  an  integral  part 
of  the  program. 

Note  that  beam-to-beam  correlation  is  not  generally 
available,  since  only  one  TL,  one  beam  pattern,  and  one  ship 
density  are  usually  run  at  a time.  It  could  be  obtained  for 
cases  in  which  beam  patterns,  ship  densities  and  TL  are  nearly 
the  same  for  the  two  steering  directions. 

2.9.8  Computer  Implementation 

SIAM  II  is  programmed  in  FORTRAN  and  is  presently 
installed  on  a NRL  PDF  11/45  computer.  It  requires  26K  core. 
Data  entry  (ships,  TL,  etc.)  is  from  files  stored  on  disk 
comprised  of  a main  file  and  separate  TL  files.  Creation  of 
these  files  can  be  tedious. 

Costs  depend  on  the  number  of  ships  per  degree, 
complexity  of  the  array  response,  number  of  subsectors,  etc. 

As  an  exantple,  25  ships.  32  copies,  2 3®-sectors  and  128  time 
points  required  10  minutes  running  time  on  the  POP  machine. 
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According  to  the  author  of  the  progra-n,  SIAM  II 
has  been  compared  with  North  Pacific  beam-noise  data  end  showed 
very  good  agreement.  Also,  see  Ref.  A-1,  I-l,  and  1-2. 


2.9.10  Significant  Advantages  and  Disadvantages 

Among  the  limitations  of  the  SIAM  II  model  are: 

(a)  The  ship  traffic  model  may  be  suitable  for 
certain  conditions,  but  does  not  allow  for 
beam-to-beam  correlation  or  the  accurate 
simulation  of  ship  tracks. 

(b)  The  array  response  and  TL  models  do  not 
account  for  vertical  arrival  structure. 

(c)  The  time  series  are  limited  to  128  sample 
points,  by  core  constraints. 

(d)  The  operation  of  the  model  is  tedious. 

(e)  The  model  is  in  general  not  structured  to 
give  the  details  of  the  noise  fluctuations 
when  complex  side-lobe  structure  or  TL 
fluctuations  are  Important. 

(f)  Only  a single  beam  pattern  is  allowed  (per 
sector)  in  one  model  run. 

(g)  The  TL  fluctuation  is  uncorrelated  in  time. 
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The  time  and  cost  considerations  are  also  Important,  but 
difficult  to  compare  with  the  other  Brute  Force  models. 

On  the  positive  side,  the  approach  Is  innovative 
and  worthy  of  further  investigation.  Some,  but  not  all, 
of  the  distinct  advantages  of  a Brute-Force  model  mentioned 
before  (e.g. , short-time  statistics,  higner-order  statistics) 
are  to  be  found  in  SIAM  II. 


Section  3 

SUMMARY  AND  RECOMMENDATIONS 

The  tables  given  at  the  end  of  this  section  sum- 
marize the  properties  of  each  of  the  nine  noise  models 
reviewed  in  this  report.  The  problem  now  is  to  choose  those 
models  or  parts  of  models  for  use  in  LRAPP  efforts:  area 
assessments,  exercise  planning,  measurement  analysis,  per- 
formance predictions  for  fleet  use,  etc.  In  lieu  of  an 
attempt  to  decide  here  what  noise  statistics  and  type  of 
model  are  needed  for  each  of  these  activities,  the  next 
several  subsections  list  some  potentially  useful  quantities 
and  identify  candidate  models  for  predicting  them.  The 
final  subsections  give  general  recommendations  for  a LRAPP 
approach  to  modeling  beam-noise  statistics  and  discuss  some 
special  problem  areas. 

3.1  ONE-DIMENSIONAL  DENSITY  FUNCTIONS  AND  MOMENTS, 

GRAND  ENSEMBLE 


Here  is  required  the  first-order  (one-dimensional) 
statistics  of  the  beam  noise  level  or  intensity,  viz.  , the 
density  function,  percentiles,  mean,  variance,  skewness. 


etc.,  but  no  temporal  statistics.  By  "grand  ensemble"  is 
meant  that  the  statistics  are  ensemble  averages  over  all 


variables  which  are  modeled  as  random,  e.g. , ship  counts, 
ship  tracks,  source  level  variations,  etc.  For  all  of 
the  Analytic  models,  such  an  ensemble  is  probably  consis- 
tent with  an  average  over  at  least  several  days  since  each 
one  of  them  defines  variations  which  would  take  that  long 
to  occur.  For  the  Brute  Force  models,  time-series  replica- 
tions covering  many  days  or  else  a number  of  random 
’'restarts'*  is  required  to  give  such  statistics. 

Hence,  the  most  efficient  approach  to  obtaining 
this  solution  is  the  USI  model.  It  is  orders  of  magnitude 
faster  than  any  of  the  other  Analytic  routines  and  does 
not  rely  on  asymptotic  limits  or  approximations  to  nor- 
mally distributed  variables.  Certain  additional  work  on 
the  model  seems  to  be  in  order,  however; 

• A computer  routine  to  automate  the  input 
(such  as  FANIN)  must  be  developed  in  order 
to  make  the  model  suitable  for  frequent 
execut ion. 

• The  Poisson  ship-count  assumption  (consis- 
tent with  WAGNER,  and  in  a sense  with  BTl. 
and  DSBN)  should  be  tested  for  sensitivity 
and  reviewed  in  light  of  data  (see  Reference 
B-2). 
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The  unusual  nature  of  the  recursive  calcu- 
lation of  the  density  function  suggests  that 
a full  documentation  of  the  underlying 
algorithms  be  reviewed  by  LRAPP. 

Moments  and  the  density  function  should  be 
calculated  both  for  noise  intensity  units 
and  for  noise  levels  (dBs).  This  should 
be  a trivial  extension  of  the  existing 
routine. 


ONE-DIMENSIONAL  DENSITY  FUNCTION  AND  MOMENTS. 
LIMITED-TIME  ENSEMBLE 


In  some  applications,  the  "grand  ensemble" 
statistics  are  not  appropriate,  and  it  is  the  short-term 
(say  12  or  24  hours)  description  of  noise  which  is  required. 
For  example,  a "grand-ensemtie"  calculation  would  not  be 
appropriate  for  simulating  the  temporal  behavior  of  measured 
noise  frcn  an  experiment  in  which  ship  positions  are  known. 
The  Analytic  models  cannot,  in  general,  produce  such  sta- 
tistics since  t.‘>ey  ensemble  over  the  entire  population  (see 
Subsection  l.l  for  further  discussion  of  this  problem). 

The  BTL  model  or  the  limiting  log-normal  approximations  of 
BTL  and  WAGNER  are  the  only  onef'»  which  produce  higher  order 
statistics  and  hence  have  the  potential  to  provide  short- 


term data. 


m 


On  the  other  hand,  the  Brute  Force  models  can 
easily  produce  i.hort-term  sample  functions  and  their  time- 
ensembled  statistics.  Any  of  the  five  models  listed  could 
provide  a set  of  density  functions  and  moments  for  short- 
term realizations. 

3.3  MINIMAL  TEMPC*RAL  STATISTICS,  GRAND  ENSEMBLE 

By  "minimal  temporal  statistics"  is  meant  the 
temporal  autocorrelation  function,  decorrelation  time,  and 
percent  of  beam-free  time.  These  can  be  important  for 
detection  analysis. 

Of  the  Analytic  models,  only  the  BTL  model  can 
provide  these  data  for  the  general  case.  Only  if  the 
Gauss-Markov  approximation  is  valid  can  the  WAGNER  Analytic 
model  produce  autocorrelation  functions  without  Monte  Carlo 
replications.  The  USI  model  can  estimate  beam-free  proba- 
bilities only,  and  the  BBN  model  cannot  at  present  predict 
autocorrelation  functions  (the  two-dimensional  density 
calculation  has  not  yet  been  implemented).  Hence,  subject 
to  an  examination  of  the  Poisson  ship-arrival  assutaption, 
the  BTL  approach  seems  to  be  the  leading  candidate. 
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HIGHER-ORDER  TEMPORAL  STATISTICS 


For  certain  applications  the  level-crossing 
properties  of  noise  (e.g. , the  probability  that  noise  re- 
mains below  a threshold  for  a given  time  period,  or  the 
waiting  time  to  cross  a level)  are  required.  These  are  in 
general  "higher-ordor"  statistics  in  the  sense  that  a two- 
dimensional  density  does  not  give  enough  information  to 
determine  them.  Among  the  approaches  available  to  predict 
level-crossing  properties,  the  Brute-Force  models  (DSBN, 
BEAHPL,  NABTAU,  SIAM  I or  II)  would  be  the  easiest  to  apply. 
In  fact,  DSBN's  current  statistical  package  described  in 
Section  2.6  calculate  level- crossing  properties.  The 
principal  drawback  is  the  expense  of  making  computer  runs 
for  multiple  replications. 

The  only  real  alternatives  to  the  Brute-Force 
approach  are  of  dubious  potential: 

• the  BTL  model  - which  involves  calculations 
equivalent  to  generating  multidimensional 
density  functions  from  characteristic  func- 
tions. or  else  use  of  the  asymptotic  limiting 
distribution  - which  still  requires  a correla- 
tion function. 


itliMMi) 


• the  Gauss-Markov  approximation  suggested 
by  WAGNER  - which  remains  to  be  verified 
for  beam-noise  applications. 

Note  also  that  the  Analytic  models  will  yield  only  grand- 
ensemble  values,  while  the  Monte-Carlo  models  can  give  short- 
term samples. 

3.5  MULTI-BEAM  AND  MULTI-ARRAY  STATISTICS 

The  prediction  of  multi-beam  or  multi-array  sta- 
tistics has  its  own  special  importance  for  current  applica- 
tions, as  well  as  its  own  special  problems.  Of  the  Analytic 
models  reviewed  here,  only  the  BTL  approach  considers  the 
general  case.  For  special  canonical  ship  scenarios.  Reference 
B-2  lists  forrolas  for  cross-correlation  functions  and  asymp- 
totic joint  densities.  The  WAGNER  model,  although  structured 
to  predict  spatial  correlation  for  omul  sensors,  has  promise 
for  multi-array  applications. 

Any  of  the  Brute-Force  models  can  be  made  to  pre- 
dict the  simultaneous  noise  outputs  from  several  beams. 

Both  DSBN  and  NABTAU  are  presently  configured  to  do  it.  and 
liSBN  has  a special  package  which  calculates  cross-correlation 
and  joint  density  functions  for  beam-to-beam  correlation 


statistics.  The  situation  is  different  for  the  case  of 
u’ultiple  array  locations.  Except  for  NABTAM,  the  Brute-Force 
models  must  be  rerun  for  each  new  receiver  location,  but 
with  the  same  ship  field.  This  is  straightforward  for  DSBN, 
BEAMPL,  and  SIAM  I,  but  not  possible  for  SIAM  II  (since  the 
ships  are  oriented  about  the  receiver  location). 

Since  the  combination  of  cases  and  model  attributes 
are  many,  Table  3-1  shows  the  models  judged  to  be  leading 
candidates  for  each  application 

3.6  GENERAL  RECOMMENDATIONS 

From  the  preceding  discussion  it  should  be 
apparent  that  no  single  model  is  the  be.st  choice  for  all 
applications.  Nor  do-^s  it  seem  tnat  pieces  and  parts  of 
the  several  models  described  cun  be  integrated  into  one 
all-purpose  routine.  However,  it  is  proposed  that  LRAPP 
take  the  following  approach; 

(A)  For  a first  look  at  grund-enstunblo . fir.st- 
order  statistics,  u.si?  the  USI  model  (after 
automation  of  input,  etc.  listed  in  Section 
3.1) 
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(B)  For  insight  about  controlling  parameters 
and  to  obtain  bounds  on  temporal  statistics, 
use  the  BTL  analytic  formulas  for  canonical 
cases.  Again,  some  kind  of  automated  system 
must  be  developed  to  generate  inputs. 

(C)  For  prediction  of  the  detailed  noise  fluctua- 
tions, for  measurement  analyses  or  detection 
studies,  use  one  of  the  Brute-Force  models. 

It  seems  that  a modest  effort  would  yield  a 
synthesized  model  employing  automated  input 
and  statistical  analysis  packages  plus  the 
best  features  of  the  five  models  listed 
(e.g.,  great -circle  sailing  from  SIAM  I, 
analysis  packages  from  DSBN,  limited-time 
and  grand  ensembles,  etc.). 

(D)  F<'r  multi-beam  correlation,  use  BTl.  formulas 
for  bounds  and  one  of  the  Brute-Forct'  models 
(of  Table  3-1,  as  appropriate)  for  details. 


3.7  ADDITIONAL  HEMARKS 

This  report  ci>ncludes  with  a brief  identification 
of  probUnns  shared  by  the  reviewed  models  and  which  require.^ 
further  study. 


• Inputs  for  Detection  Simulation 

None  v>f  the  models  reviewed  can  be  used  to 
simulate  the  real  deiecti\»n  prtH'oss,  i.e,. 
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to  predict  signal-plus-noise  in  one  frequency 
bin  and  an  "estimate"  of  noise  (e.g.,  value 
in  another  bin  or  average  over  other  bins,  etc.). 

Noise  Statistics  for  Time-Dependent  Array 
Responses 

In  actual  applications,  the  array  may  well 
be  towed,  so  that  then  its  location  changes 
in  time  and  its  response  changes  in  time 
(distortion).  None  of  the  models  deal  with 
this  case.  The  beam- free  times  and  other 
noise  fluctuation  properties  are  expected 
to  be  affected  significantly  by  this  array 
motion. 

Ship  Infoimmtion 

The  most  difficult  and  time-consuming  task 
associated  with  producing  beam-noise  statis- 
tics seems  ti)  be  the  preparation  of  surface- 
ship  densities  and,  t'specially,  .ship  lanes. 
Although  some  automated  routines  exist  (e.g.. 
KANIN,  the  input  for  FANM , ut  NORDA)  they 
are  not  widely  u.sed,  and  further  development 
and  eviiluatlon  are  indicated. 

In  addition  to  the  ship- local  ion  problem, 
there  remains  a significant  uncertainty 
in  the  prediction  of  .source  levels  and 
d I rec  1 1 V 1 1 y . Uecen  I nu'astiremen  t s st»ou  I d 
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the  radiated  noise.  For  the  present,  how- 
ever, the  Ross  and  Alvarez  spectra  of  Ref, 

E-2  and  the  BBN  statistical  model  of  Ref. 

C-3  (which  accounts  for  discrete  components 
of  the  spectrum)  are  the  best  data  available. 

• Noise  Model  Evaluation 


If  convincing  validations  (with  measurements) 
of  mean-noise  models  are  scarce,  then  those 
for  beam-noise  statistics  models  are  very 
rare  indeed. 

• Wind-Dependent  Noise 

Little  attention  has  been  given  here  to  the 
modeling  of  the  statistics  of  wind-related 
("surface"  or  "sea-state"  or  "wind/wave") 
noise,  either  omni  or  at  the  array  output. 
Such  noise  can  be  important,  oven  at  low 
frequencies,  and  the  problem  warrants  fur- 
ther investigation. 
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Summary;  ANALYTIC  MODELS  (I) 


Approach 


The  method  involves  the  direct  calcu- 
lation of  one-  and  two-dimensional  pro- 
bability density  functions  by  Fourier 
inversion  of  the  characteristic  func- 
tions. The  important  calculation  is 
of  the  characteristic  function  when 
shipping  traffic  and  source  levels  are 
random  variables.  The  model  is  struc- 
tured to  provide  "main  beam  noise"' 
only  - i.e.,  noise  from  a fi.xod  azi- 
muthal sector. 


Ship  Source  Intenslt 


Ship  source  level  is  a random  variable, 
independent  from  ship  to  ship,  selec- 
table by  class,  but  always  independent 
of  time  and  azimuth  and  vertical  angle. 


The  noise  process  is  treated  as  "shot- 
noise  process'"  by  the  assumption  that 
ships  travel  on  paths  and  have  Poisson 
arrival  times.  Calculation  of  multi- 
dimensional characteristic  functions 
and  moments  use  shot-noise  assumption. 
Fourier  inversion  or  use  of  log-normal 
limiting  distribution  Is  often  needed 
to  obtain  statistics. 


Same  as  BBN.  but  all  ships  on  a path 
have  the  same  distribution  of  source 
intensity. 


The  key  to  the  approach  Is  that  the 
number  of  ships  in  a geographic  domain 
la  a Poisson  variable  I'apatlal  Pois- 
son process" 1.  Then  the  noise  is  a 
single-time  sample  frisn  "generalized 
Pol.S8on  process."  The  one-dimensional 
density  is  calculated  with  a special 
Iterative  algorithm  for  which  geo- 
graphic domains  are  grouped  accord- 
ing to  contribution  iper  ship)  to  noise 
intensity.  The  method  avoids  charac- 
teristic functions  and  direct  convolu- 
t tons. 


Same  as  BBN,  but  all  ships  counted  by  the 
Poisson  variable  have  the  same  distribu- 
tion of  source  Intensity. 


ffAONKH 


Two  different  formulations  flnmnded 
process”  and  spatial  1‘olsson  distri- 
bution of  ships)  lead  to  a "general ized- 
Poisson  process"  characteristic  func- 
tion. Also,  the  approach  Includes  a 
brutt-force  method  for  obtaining  sam- 
ple paths.  A central  limit  theorem  and 
Monte  Carlo  simulation  to  obtain  corre- 
lation functions  leads  to  conclusion 
that  omni  noise  level  idB's)  la  appros- 
Imately  a Qausa -Markov  proce.sa.  A 
formula  for  the  apatial  correlation 
of  noise  la  obtained  for  the  general 
Poisson  proceaa.  The  model  is  struc- 
tured to  produce  <utini  noise  statistics 
at  distributed  locations. 


Same  as  5DN.  but  main  reHUlt.H  asnume  all 
ships  have  the  same  source  lutenslty 
distribut ion. 


1?.-- 


Sunaiary:  ANALYTIC  MODELS  (11) 


Model 

Ship  Locations  and  Motion 

Sound  Transmission 

1 

BBN 

1 

1 

Ships  are  assumed  to  travel  in  lanes  or 
routes.  The  number  of  ships  in  a lane 
of  a certain  class,  the  source  level, 
initial  position,  and  course  are  random 
variables.  The  speed  is  a constant  for 
the  class,  while  source  level  distri- 
bution is  the  same  within  a class.  Ships 
travel  at  constant  course  and  speed  once 
initialized;  "Rectangular  Sailing." 

Transml.'s^ ton  loss  is  deterministic  and 
defined  according  to  azimuthal  sector 
and  ran >>8.  There  are  no  restrictions 
on  the  amount  of  detail  permitted.  The 
effects  of  vertical  arrival  structure, 
coherence,  etc.,  are  not  Included. 

1 

BTL 

Ships  travel  on  "paths,"  which  may  be 
grouped  in  lanes  or  as  isotropic  fields. 

On  a path,  ships  arrive  at  the  main  beam 
steering  angle  according  to  a Poisson 
rule.  Furthermore  ships  on  a path  have 
course,  speed,  and  source  intensity 
drawn  from  the  same  distribution.  Once 
initialized,  a ship  travels  at  constant 
course  and  speed; 

"Rectangular  Sailing." 

i 

[ 

i 

i 

Same  as  BBN 

1 

1 

. 

j 

OSl 

1 

Since  the  model  does  not  produce  tomporsl 
statistics,  the  ships  do  not  move  in 
time.  Input  is  the  mean  number  of  ships 
(by  class,  if  desired)  in  -‘acn  of  many 
geographic  domains.  The  number  of  ships 
in  each  domain  and  class  la  then  treated 
as  a Poisson  variable,  and  those  ships 
are  assumed  to  have  the  same  source 
level  distribution. 

Same  as  BBN,  except  that  a TL  fluctuation 
distribution  can  be  defined  at  input  for 
each  geographic  region  in  which  ships 
are  defined.  This  TL  fluctuation  variable 
l»  combined  with  the  source  level  variable 
in  the  model  calculations.  The  fluctuation 
value  is  independent  from  ship  to  ship. 

j 

1 

VAQNEB 

The  salps  ar«  distributed  initially  ac- 
cording to  a spatial  Poisson  process, 
l.e. . the  number  of  ships  in  a region 
is  a Poisson  variable  whose  mean  i.y  pro- 
portional to  the  area  of  the  region,  or 
elte  aa  a hounded  procesa  in  which  a 
fixed  number  of  ships  are  distributed 
uniformly  in  a bounded  region.  In  each 
case,  the  course. speed  and  aource  level 
are  random  variables,  choaea  independently 
for  each  ship  from  tbs  same  distributions; 
once  selected  these  variables  remsln  con- 
stant in  tine.  For  the  Poisson  process, 
ships  are  defined  on  an  unbounded  region, 
tor  the  bounded  noise  process,  ships 
rsflect  from  boundaries  of  the  douin. 
Sailing  la  "rectangular." 

Sane  as  BBN.  A single,  azlmuihai 1 y- 
independent  TL  is  used  at  present 

Model 


Sumury:  ANALYTIC  MODELS  (III) 


Receiver 


The  receiver  Is  a horizontal  array  with 
fixed  location  and  deterministic,  time- 
independent  response.  All  energy  is  as- 
sumed to  arrive  as  plane  waves  in  the 
horizontal  plane  - perfectly  coherent 
and  undistorted.  The  model  is  struc- 
tured to  provide  output  (or  the  main  beam 
only  U-e.,  noise  from  an  azimuthal  sec- 
tor). Contrtbut ion.s  from  di  (ferent  sources 
are  added  incoherently  (random  phase). 
Extension  to  predict  sldelobe  levels  is 
possible  - but  computationally  expensive. 
The  array  is  assumed  to  be  perfectly  stralgl 
horizontal  and  ideal  in  its  response. 


Output 


Fully  ensembled  one-dimensional  density 
and  moments  of  noise  intensity  in  one 
azimuthal  sector.  Eventually  the  model 
will  produce  two-dimensional  densities 
(in  time),  autocorrelation  functions, 
etc . 


I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 


BTL 


Same  as  BBtf.  except  that  noise  from  the 
sidelobes  can  be  predicted.  The  model 
has  also  been  structured  to  produce  multi- 
beam, multi-array  noise  statistics. 


Fully-ensembled  k-dimensional  densities 
and  moments  are  available.  Computations 
become  difficult  for  k > 2.  For  special 
cases,  there  are  formulas  for  momenta 
and  correlation.  The  model  can  produce 
beam  to  beam  or  multiarray  joint  statis- 
tics. If  a certain  asymptotic  limit  applies, 
the  noise  intensity  Is  modeled  as  a multi- 
variate log-normal  variable.  In  that  case, 
the  calculation  of  higher  order  statistics 
is  greatly  simplified. 


US  I 


Same  as  BTL.  However,  a random  array - 
response  fluctuation  for  each  ship  group 
can  be  imposed.  It  is  treated  in  the  same 
way  as  the  source  level  and  TL  fluctuations. 


The  fully-ensembled  one-dimensional  density 
function  (with  1.5  or  3dB  resolution)  and 
moments  of  main-beam  and  side- lobe  noise 
Intensity  are  predicted.  Model  provides 
Infomation  on  which  ships  are  tmfKirtunt 
to  mean  and  variance.  No  higher  order 
statistics,  except  percentage  of  beam 
free  times,  are  available. 


WACNER 


The  receiver  la  an  >)«nl  phono  lor  phones' 
with  fixed  location  and  deterministic,  time- 
independent  response.  Extension  to  the 
horizontal  army  case  is  possible.  The 
model  Is  structured  to  produce  multi- 
sensor  noise  statistics. 


The  brute-. jfce  implementation  allows  for 
the  Monto-Oarlo  calculation  of  higher 
order  statistics.  Otherwise,  the  aaidel 
produces  the  characteristic  function,  one- 
dimensional  density  and  moments,  fully- 
ensembled.  No  analytic  formulas  for  higher 
order  (tinel  siatisiica  are  given.  The 
spatial  correlation  function  is  written  in 
integral  form  and  approximations  are  li.sted. 
j Ihen  the  Ciiusn-Uarkov  approximation  i»  valid, 
higher  ordiM-  statlstica.  level -crosaing  re- 
sults and  time-series  simulations  are  readily 
available. 


_ijnij_iji,mtTtfiig  iiTTi ^ 
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Sumuiftry  ANALYTIC  MODELS  (IV) 


Oc^nput Ational  Elfuueiiey 

• AU  tnpvit  la  manual. 

• Computer  oaloulatton  ot  oharai'torlst lo  functions  and  Fourier 
inversion  (or  one-d)mensloiial  density  (50-60  dB  dynamic  ranee) 
takes  about  one  CDC-6400  minute  i$lO). 

• Two  dimensional  statistics  are  expected  to  be  calculated  at 
mucU  greater  expense. 


• Same  as  BBN  except  that  the  form  of  the  characterist Ic  (unc- 
tion and  a .smaller  dynamic  range  (30  dB)  yields  less  e.xpen- 
sivo  routine. 

• BTL  IS  structured  to  calculate  higher  order  and  multiarray 
stat 1st ics. 

• When  the  asymptotic  limit  applies,  a very  efficient  model 
results. 

« Typical  calculation  ot  one -dimensional  density  functu'ii  on 
L'NIVAC  1106  takes  dOK  core  and  6 minutes,  per  shipping  lane 


• Alt  input  IS  manual,  although  a FANlN-type  routine  would  help, 
since  only  ship  den.-uties  are  used. 

• TL.  SI..  AO  ,ind  .ship  domainti  are  grouped  manually. 

• Actual  computation  of  density  function  costs 
less  than  $1.00. 


« All  input  IS  manual.  Ship  input  i.s  not  difficult  since  Initial 
! distribution  lu  uniform  and  course*  and  speeds  are  given  by 

single  distributions 

! s Tbs  analytic  s^’l^tlon  for  one-diaensionsl  density  requires 

Fourier  inversiv*n  of  chsracirri st  iv*  functions  illke  BBN  or  BTL). 

• I £ (lauss-Marko*.  appro vinat  too  is  lalivl  Wagner  has  mam  closed 
i fora  result*  and  computer  simulativ'n  iivil*. 

I 
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SuiHury:  BRUTE  FORCE  MODELS  (I) 


Model 

Ship  Source  Intensity 

Ship  Locations  and  Motion 

^ BEAUPL 

Ship  source  levels  sre  deterministic 
constants  for  each  class  of  ship  and  are 
Independent  of  angle  and  time. 

Same  as  DSBN,  except  that  all  ships  In  a 
lane  have  same  deterministic  constant 
speed  and  source  level. 

OSBN 

Ship  source  level  is  a random  variable 
depending  on  class,  speed  and  length. 

The  formula  is 

SL  ■ K ♦ 60  log  (Speed)  > 20  log 
(Length),  where  K Is  a constant  for  each 
class.  Speed  and  Length  are  uniform 
variables  drawn  Independently  for  each 
ship.  The  levels  are  constant  in  time 
and  independent  of  angle.  Extension  to 
Include  time  or  angle-dependent  level  Is 
possible,  but  adds  to  computation  load. 

Ships  travel  on  paths  with  Poisson  arrival 
rule  for  Input  ship  path  density  (like 

BTL).  Route  envelopes  on  lanes  are  defined, 
and  within  them  ship  properties  (course, 
speed,  level)  are  drawn  from  the  same  dis- 
tributions, independently.  Once  initialized, 
ships  sail  (rectangularly)  at  constant 
course  and  speed.  Lanes  are  initialized 
on  circular  arcs  at  edge  of  basin. 

NABTAM 

1 

All  surface  ships  have  the  same  constant 
source  level,  derived  from  Ref.  E-2,  de- 
pending only  on  frequency. 

1 

Initial  ship  locations  are  a deterministic 
input.  Initial  courses  and  speeds  can  be 
user-specified  or  drawn  from  uniform  and  nor- 
mal distribution,  respectively. 

Once  initialized,  ships  travel  on  tracks 
with  dLAT/dt  and  dLON/dt  constant.  If  a 
ship  leaves  the  basin,  it  is  not  replaced. 

SIAM  I 

Source  Intensity  is  a random  variable 
tor  each  ship,  determined  from  distribu- 
tion for  ship  class,  and  depending  on 
frequency,  and  uniform  or  normal  vari- 
ables for  speed  and  length.  Default  val- 
uee  are  the  same  as  for  OSBR. 

Source  inteneity  does  not  depend  on  time 
or  tranaalselon  angle. 

j 

Ship  locatione.  courses  and  speeds  are 
initialized  by  class  either  determinis- 
tically (to  match  experimental  data) 
or  randomly  (usually  uni  fora). 

Ooce  determined,  the  ships  travel  on  great 
circle  paths  at  constant  speed  and  with 
constant  source  level.  When  a abip  reaches 
ths  basin  boundary  it  Is  reflected  speculsrly 
(as  la  the  VAGRER  bounded  cotse  process). 

As  asay  as  250  ships  are  allowed 

SIAM  tt 

1 

1 

Same  as  SIAM  I 

Initial  ship  courses  and  spweds  are  chosen 
from  near-normal  distribution,  but  then  re- 
main constant  in  tine.  Initial  ship  loca- 
tions are  uniform  in  subsectors. 

Ships  are  tracked  in  time  on  straight-! tne 
paths  with  reflections  from  subsector  boun- 
daries (see  2.9.3  for  details) 

1 ’ 

Sunoary;  BRUTE  FORCE  MODELS  (II) 


BEAMPL 


Sound  Transmission _____ 

Transmission  loss  is  a deterministic  in- 
put and  is  assumed  to  be  independent  of 
azimuth  and  time.  There  are  no  restric- 
tions on  the  amount  of  detail  permitted. 
The  effects  of  vertical  arrival  struc- 
ture. coherence,  etc.  are  not  included. 


Receiver 

The  receiver  is  a perfect  horizontal  array 
with  fixed  location  and  deterministic, 
time-independent  response.  In  fact,  the 
response  in  the  horizontal  is  simply  unity 
on  the  "main  beam"  and  zero  elsewhere.  All 
energy  is  assumed  to  arrive  as  horizontal 
plane  waves  - perfectly  coherent  and  undis- 
torted. 


■?! 


.VABTAl! 


SIAM  I 


Same  as  BEAMPL,  except  that  TL  can  be 
given  as  a function  of  azimuth  and  ver- 
tical arrival  angle. 


Transmission  loss  is  calculated  with- 
in the  model,  using  a geometric- 
acoustics  approach  without  diffrac- 
tion or  coherence  corrections,  for 
a range- Independent  environment. 

The  TL  and  vertical  arrival  struc- 
ture are  calculated  for  each  of 
several  receiver  locations  to  all 
source  ranges.  The  TL  is  deter- 
ministic and  time-invariant. 


A deterministic,  time- independent 
TL  as  a (unction  of  range,  azi- 
muth. frequency  and  receiver  depth 
is  required  for  input.  Up  to  30 
frequencies.  39  depths,  and  10  azi- 
muths are  .allowed.  At  each  time 
sample  an  independent  draw  from  a 
fluctuation  distribution  is  added 
to  the  TL  for  each  ship  source. 
Also,  there  is  a special  'logarith- 
mic" routine  for  interpolating  TL 
in  range  and  azimuth,  and  all  geo- 
eaetry  is  spherical. 

The  program  can  generate  TL  of 
form  .A*BlogR  internal  ly  - as  an 
option . 


The  model  assumes  an  array  at  a fixed  loca- 
tion with  beam  pattern.  Subroutines  generate 
canonical  patterns  for  shaded  and  unshaded 
line  arrays.  All  energy  Is  assumed  to  arrive 
as  plane  waves,  distributed  in  vertical  an- 
gle. but  perfectly  coherent  and  undistorted. 
The  model  can  account  fully  for  the  vertical 
multipath  arrivals  and  their  effect  on  array 
response  (or  angles  off  broadside 


Receiver  locations,  depths  and  beam-patterns 
are  deterministic  Inputs.  Internal  routines 
can  calculate  beam  natterns  for  standard 
line  array.s.  The  three-dimensional  response 
(including  vertical  multipath  arrivals)  is 
calculated  as  in  DSBN 

.Also  there  is  a version  of  the  model  which 
performs  "near-field"  corrections,  l.e., 
it  calculates  array  responses  for  curved 
wavefronts  from  nearby  sources. 

Same  as  BEAMPL  except  that  multiple  beam 
patterns  can  be  treated  in  one  program 
execut Ion . 


I 

1 1 
I I 


SIAM  n 


The  model  requires  as  input  a 
table  of  TL  versus  range  (or  as 
many  as  Id  azimuths,  one  fre- 
quency. and  one  receiver  depth. 

A fluctuation  component  and 
A*BlogR  routine  are  available, 
as  10  SIAM  I. 


Receiver  locatu^n  and  depth  are  fixed  .1 
single  beam  pattern,  depending  onSv  on 
azimuthal  arrival  angle,  is  an  input  and 
defines  the  t ime- independent  array  re- 
sponse. 


tHy  I*' 


Sujnn«.ry  ; BRUTE  FORCE  ROOELS  (III) 


Model 


Output  Antlvsls 


Basic  Output 

The  basic  model  output  is  a sampled  time 
series  of  beam  noise  for  extended  tine 
periods  (hours  or  days).  Monta-Carlo 
simulation  requires  that  a number  > f such 
replications  be  calculated. 


The  analysis  package  generates  the  statis- 
tics for  the  time  series.  Including: 

(a)  Time-averaged  series  (intensity  average) 
for  each  replica. 


I 

i 

I 

I 

I 

I 

I 


BEAMPL 


The  model  also  produces  a complete  his- 
tory of  each  ship. 


(b)  Histograms,  moments,  percentiles  for  each 
replica. 

(c)  Time-averaged  autocorrelation  function 
and  its  Fourier  transform  (dB's)  for 
each  replica, 

(d)  Plots  of  (a),  (b).  (c). 

(e)  A stationarlty  tost  Is  performed  by  find- 
ing moments  and  percentiles  for  subinter- 
vals of  the  time  series. 


(f)  Level  crossing  statistics  - percentage 
of  time  noise  is  above  a given  value. 


USBM 


Same  as  BEAMPL.  except  the  series  repre- 
sents full  array  response  - both  from 
main  beam  and  sldelobes.  The  Bx>del  can 
also  calculate  the  responses  for  several 
beampatterns  at  once  (l.e..  duo  to  the 
same  ship  field),  so  that  beam-to-beam 
correlation  studies  can  be  made.  The 
model  gives  ship  histories  plus  list  of 
ships  on  main  beam  and  on  sldelobes. 


Same  as  BEAMPL  {(a)  through  (f)1.  plus: 

(g)  Two-and  three  dimensional  histograms 
(in  time  or  beam),  plus  moments,  cross- 
correlations.  etc. 

(h)  Comparison  of  Intensity  distribution 
«lth  best  fit  chi-square  (central  and 
non-central),  log-normal,  etc.  distri- 
butions. Kolmogorov  test  results  and 
plots  are  displayed. 


(1)  Ensemble  density  and  ^statistics  over 
all  replications  at  one  time  period 
are  found. 


(J)  E.xtensive  level-crossing  reault.s.  In- 
cluding distribution  of  time  intervals 
for  ehlch  noise  exceeds  a level  and 
waiting  times  to  cross  a threshold. 


(k)  Cotaplete  history  of  beam-free  and  beam- 
occupied  times. 


1 

I 

I 

I 

I 
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SABTAM 

Basic  output  consists  of  noise  time 
series  for  each  receiver  location  and 
bean  pattern  at  one  frequency.  Monte- 
Carlo  replications  are  not  generated. 

Mo  analysis  package  Is  used  at  present,  but 
a standard  titae-series  statistics  package 
could  be  applied. 

SIAM  I 

Basic  output  Is  noise  titae  series  for 
each  beae-pattern.  frequency,  and  depth 
for  extended  time  periods  and  multiple 
replications.  Ship  location  histories 
are  also  available. 

The  analysis  package  provides  the  onc-4t*ea- 
oional  density  function  and  nocwnta  for  the 
ensemble  of  rcpl Icalioni . Although  autvicorre- 
i at  Ion  functions  can  be  found,  their  interpre- 
tation depends  on  the  TL  fluctuation  model. 
Beamfree  and  occupied  percentages  and  time 
histories  are  outputs. 

SIAM  11 

Basic  output  consists  of  32  realisations 
of  beam-ncise  time  series,  for  the  single 
beam-pattern,  frequency,  and  receiver 
location,  and  up  to  129  time  steps  each. 

An  analyais  package  compute*  noise  histograms, 
moments,  and  the  median. 

Costs  depend  on  nuir.  ai  ships,  nuaber 
of  beam  ps.tterns.  replications,  time 
period,  amount  of  analysis,  etc.  Examples 
are  (CDC  6400  computer): 

• 200  ships,  12  hours.  9 beam  patterns, 
no  analysis:  $9/ replication. 

e Analysis  of  10  10-hour  replications 
including  (a),  (b).  (c).  (d),  (e). 

<f)  (h).  (i>.  (J).  (k):  $20. 

Core  requirement  for  CDC  6400  Is  less 
than  40K  words. 


} • 


.lADTAy 


Program  is  tverlaid  (cr  ainiccaputer  appli- 
cation. and  cor*  is  lOK  words  iif  not  over- 
laid. t3K). 

E.Tample  of  cost  one  tis^  step,  large  num- 
ber of  ships,  five  beam  patterns  $0.T$. 
on  CDC  6400.  Note,  cost  includes  TL  calcu- 
lation. 


A typical  run  on  the  CDC  WOO  requires  5-10 
minutes  computer  time. 


SIAM  U j 


Costs  depend  on  ships,  array  response,  sub 
sectors,  etc.  fot  exarqile.  2.5  ship.s.  02 
realisstScas.  2 sectors,  125  time  points 
per  realisation  requires  10  minutes  on  the 
MM*  11/45  Polluter 


A-2. 


A-4. 


A-6. 


B-2. 


B-3. 


C-l. 


C~2. 
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